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Zum Problem des ebenen Spannungszustandes im kreis- 
zylindrischen spröden Rohr unter konstantem 
Innen- und Außendruck 
Von E. Glatzel in Marl-Hüls und H. Schlechtweg in Kettwig/Ruhr 


Unter Benutzung des nichtlinearen Verformungsgesetzes spröder Stoffe wie Gußeisen, Sandstein usw. 
wird im ersten Kapitel die Differentialgleichung 2. Ordnung für die Radialspannung 0, näherungsweise mit 


ungefährer Fehlerabschätzung gelöst für Fälle, in denen m < —1 bzw. - >2 ist. Im zweiten Ka- 


r 2 
pitel handelt es sich um eine andere Lösungsmethode, die sich einer Differentialgleichung für die Radial- 
verschiebung bedient, zu deren Aufstellung das benutzte Verformungsgesetz näherungsweise nach dem 
Spannungstensor aufgelöst wird. 


The paper is based on the non-linear law.of elasticity of brittle materials, such as cast iron, sandstone elc. 


In the first chapter, the differential equation of the radial tension o, is solved approzimalely, for = <—] 
E 


[0] 
and nn > 2, and an approximate estimation of the error is given. In the second chapter, another way of 


f 
solution is treaied, that uses the differential equation of Ihe radial deviation. For this purpose, the law of 
elasticity is approximately solved for the tensor of tensions. 


A lapplication de la loi d’ölasticite non lineaire de maltieres cassantes, comme fer de fonie, gres etc., au 

premier chapitre l’Equation differentielle de second ordre pour la tension radiale o, est resolue approzimalive- 
< . . [27 O9 

ment avec Evaluation approwimalive des erreurs pour les cas, dans lesquels =. <= 1 0u = >= 


r r 
Au second chapitre il s’agit d’une autre methode de solution, qui se sert d’une Equation differentielle pour la 
deviation radiale, a la formation de laquelle la loi d’elasticite en question est r&solue approwimativemen! 
second le tenseur de tensions. 


B nepsoä ruaBe npenuaraeMoü FaboTk IIPHÖNKKeHHO pemaeTca AußpepeHumanbHoe ypaB- 
HeHHEe BTOPOTO IOPAAKa AA PANHaJIBHOTO HAIPAKEHHA C, C IIPHÖNHSHTEIIBHOH ONEHKOH OIMIHÖOR 


[07 [0] u “ 
Ha, CIIYYaH, KOTA& zn <—1l um > >2; AA 9TOTO OPUMeHAeTcA HelIHHeÜHRH 3aKOH Nebop- 
r r 


MaUuH XPYIIKUX MaTepuaoB, KAK YyTyH, NecyaHuk u T.ı1. Bropan r1aBa NocBAmeHa APpy- 
TOMy MeTony pemeHusn, UCHONB3YIOIeMy OAHO AubhepeHnmanbHO® YP&BHeHHe IA PaAAHaJIb- 
HOTO CABHuTaA. Ü 3Toi MelibIo BEIIIeYIIOMAHYTEIÜ 3AKOH MedopMalun IIPHÖNNKEHHO PemmaeTcH 
OTHOCHTEABHO TEH30PA HAIPAKEHHA. 


IaKR aBIitel 
Näherungslösung einer Differentialgleichung für die radiale Spannung 
1. Das nicht-lIineare Verformungsgesetz 
Das H o o ke sche Gesetz, wonach die Deformationen proportional den Beanspruchungen 
sind, gilt für gewisse spröde Körper, wie Gußeisen, Sandstein, Marmor u. dgl., nicht mehr. 
Sie folgen eher dem folgenden nicht-linearen Verformungsgesetz !), ?), 
sie eg; O-4t ee 
1 — BIP — OYP + ZUR) 
wobei nur die positiven Wurzeln zulässig sind. 
Es bedeuten: 


u die Querkontraktionsziffer 
E den Elastizitätsmodul 
B, CO und ZL weitere Materialkonstante 


_*) Die Arbeit ist als Dissertation der Philosophischen Fakultät zu Köln von E. Glatzel entstanden; 
Kennziffer D 38. 
1) Schlechtweg, Über ein allgemeines Elastizitätsgesetz spröder Körper. Z. angew. Math. Mech. 14 
1934), S.1. 
2) Schlechtweg, Das Verhalten spröder Körper unter hohen Drucken und der spröde Bruch. Ing.- 
Archiy 4 (1933), S. 263. 
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& den Deformationstensor 
|S| die Summe der drei Diagonalkomponenten von © 
& den Deviator: © =& —4+|6|E&; 
€ den Einheitstensor 
den Spannungstensor 
|P] die Summe der drei Diagonalkomponenten von ® 
’ den Deviator: !=B-+|PiE; 
Indem man die Quadrate seiner Komponenten addiert, entsteht das skalare Produkt des 
Tensors ®’ mit sich selbst: (P’)*. 
Zur Abkürzung sei gesetzt: 


I=- BP +OYBP-IVRPIR| - -- 2... 0% 
J kann man?) die ,‚Anstrengung‘‘ des Materials bezeichnen ; J ist eine Invariante des Spannungs- 
zustandes. 
Für den Tensor © liefert (1) die Beziehung: 
— u 5) 1 1— 2 Ne: 3): 
© = E B ah er E #677 ee)” 


Ist bei einem linearen Zug-Druckversuch von den Spannungen nur o, + 0, so folgt für 
die Dehnung in gleicher Richtung: 


ve 
ae 


ET | 
re ke ne 


1-(F84 oy3)Ie- = EV alone | 


Das obere Vorzeichen gilt im Druckversuch: o,< 0; und das untere im Zugversuch: 0, >0. 


(Plastischer 
Verschiebungsbruch,) 


a 


Bild1. Verlauf der Anstrengung J beim linearen Zug-Druckversuch nach Gleichung (4). 


er en. a :, 
(B =4.10-kg-iemt; 0 y: =1:10kgketem; L 3 —1,5-.10-kg ? eme.) 


Nun tritt bei Zug- und Torsionsversuchen ?) der „‚spröde Bruch‘ ein, wenn die Anstren- 


gung J einen Wert zwischen 0,62 und 0,72 erreicht, hierdurch kann wegen Ei u ee in (4) 


4 2 
1—2 
stets 3 2 gegenüber 1 mit einem Fehler von höchstens 12%, als innerhalb der Fehler- 
grenze bei spröden Zugversuchen, vernachlässigt werden. — Bei Druckversuchen nimmt J 


negative Werte an, die bedeutend größer sein können. In Abb. 1 ist der Verlauf von I nach (4) 
für ein bestimmtes Wertetripel (B, €, L) skizziert. Der bei hohen Drucken auftretende „pla- 


stische Verschiebungsbruch‘“, der wohl von Y($’)? abhängt, tritt erfahrungsgemäß ein, bevor I 


°®) Schlechtweg, Der Spannungszustand in der rotierenden spröden Scheibe. Ing.-Archiv 5 (1934), S. 21. 
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den Wert I,.;, erreicht hat. / kann dann — 3 und weniger betragen. Der prozentuale Fehler 


rl: sche 1— 
bei einer Vernachlässigung von -—_ I gegen 1 kann also bei hohen Drucken und u+ 
bedeutend größer als 12% sein. 
Trotzdem gibt: 
Gt 
er ee (2), 
neben dem Zug- auch den Druckversuch richtig wieder 2). Die Forderung aber, daß ! — = 


in (4) praktisch keine Rolle spielt, ist SIASHDEITDIENT mit der Schreibweise des nicht-linearen 
Verformungsgesetzes in der Form: 


a ey 
Se ee Te 
denn das liefert für © die Tensorgleichung: 
| = [u ) 1 
SS — E Bu ee, . . . . . . . . . . (2); 


welche ohne Vernachlässigung auf (5) führt. 


In (6) hat u wieder — im Gegensatz zu (1) — die Bedeutung des absoluten Verhältnisses 
von Querkontraktion zur Dehnung. 

Der Bequemlichkeit der Rechnung halber wird im folgenden (6) bzw. (7) als das für 
spröde Körper gültige Verformungsgesetz ee L und ‚‚nicht-lineares Verformungsgesetz‘“ 
genannt. 

Bei Zugversuchen, bei denen J stets positiv ist, kann L = 0 gesetzt werden !): ‚Die 
Zugfestigkeit ist so gering, daß das mit L multiplizierte Glied nicht zur Wirkung gelangt und 
vernachlässigt werden kann.“ 

Verallgemeinert man nun die bei Zug- und Torsionsversuchen gefundene Bruchbedingung, 
indem man ihre Gültigkeit auch für alle mehrachsigen Spannungszustände annimmt: 


Der „spröde BruchHt ri # betlernem Werft von Jzwischen 
0,62 und 0,72 ein; 
so läßt sich innerhalb der Fehlergrenzen bei Untersuchungen an spröden Körpern L=0 
setzen im Falle |®| > 0 mit der für praktische Werte von B, C', L erfüllten Voraussetzungen: 


BY 


a IB Pe. 
MR sera a Ben ua A N = 


denn dann erfolgt bei genügend hohen Spannungen der ‚‚spröde Bruch‘. 


3 
Damit 7 in (4) überhaupt negativ wird, muß mindestens gelten: B> (C V3 . Auf Grund 
bisher bekannter Werte kann man sogar B> (Ü voraussetzen. x 


2, Die Differentialgleichung für die-radialerSpannung. 

Wird ein Rohr, dessen Länge groß ist gegenüber seinem Durchmesser, mit radialem 
Außen- und Innendruck 9, bzw. p; beansprucht, so liegt in guter Näherung ein ebenes Span- 
nungsproblem vor. 

R, und R, seien Außen- bzw. Innenradius. 

In Zylinderkoordinaten 7, 9, z lauten bei verschwindenden Schubspannungen die Gleich- 
gewichtsbedingungen ®): 


de, %— 0% 
A EN AR ERSTE NEN TEN, 
Fa (9) 
Für die axiale Spannung ist dabei vorausgesetzt: o, =. ; 
In dieser Arbeit möge nur der Bereich interessieren, in dem für die radiale Spannunggilt: 
et. 


Die Komponenten des Verzerrungstensors e, und e, sind durch die radiale Verschiebung s 
‚miteinander verknüpft *): 


4) Vgl. z.B. Madelung: Die math. Hilfsmittel d. Physikers. Berlin 1925; auch: Handbuch der 
Physik, Bd. VI, Berlin 1928, (Artikel Trefftz, S. 81). 
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ferentialgleichung für s (siehe 11. 
Weg zu einer Lösung känn entweder über eine Dif 

Sr benötigt man die Auflösung von (7) nach $, oder über eine Differentialgleichung 
für die radiale Spannung führen. Zur letzten gehören lineare Randbedingungen: 


für =, alt ae 
ür v=R, gilt: :, = — Pa | 
Im folgenden seien = r 

ee do yes 

En, und % a; | 
unabhängige bzw. abhängige Veränderliche. Sie sind ebenso wie die Abkürzung (ey u 
| rt . (13) 

p —= 
Ruo— Ri 


dimensionslos. yund o, werden jetzt als Funktionen von & betrachtet. Mit Hilfe von (9) wird :°) 


2 
2=,+,=(p+R)d+2o;, Y@F=Vz=:| 


. (14) 
= +3, = VRR 3Plo, +30 
und nach (7): : : 
lt I 
Be ee Bar; | 
ER al +2), 0, #42) +20))-7 i 14 5 (18). 
Pt 4 a IHR 3 
1l+u ; 1 
= (os 6: (p+2)o,+ a r 
Die letzte Gleichung liefert: 
2= let) EL) 
Eliminiert man aus. den beiden anderen ne y, so folgt: 
J ’ 
BR -a+I;| an 


14.1 BB —y@H Al]: 


bei Verwendung der Abkürzungen: 


B=B:; = 0V5; ee 


Für einen elastischen Körper gilt die mit dem Index 7 gekennzeichnete ‚Hooke sche 
Lösung 


Pa=0 PBa=— 28; 4 | 
F 2 2 9 f 119). 
- = +2); = ga) 
0,H Se une OyH 1 ?\p+x k 
Die Randbedingungen bestimmen die Integrationskonstanten. 
Im Falle 9, = 9; = Pi. folgt aus (17): 
0, — 09 — — Pia: | 
pn — l—u)Pia(P+ 2) . (20). 


E{l tr (2 ß et y) Pia 24 PraV Pia} | 
Die Lösung ist unabhängig von R,und R; und gilt somit auch für den Yolzyia 
Der Bruch ist dann ein plastischer. 


Die Tatsache, daß bei gleichem Außen- und Innendruck S —=1 gilt, ist auch in elasti- 


- 
schen Körpern erfüllt. Aus der Ho oke schen Lösung folgt weiter: 


mit 9, >pmr elt: - Bi 
# (21). 
mit m <p gilt: en a! 


5) Die Ableitung nach & wird durch einen Strich über der abgeleiteten Größe gekennzeichnet. 


} 


er 


2 ae rt 
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Es entspricht dem physikalischen Empfinden, die Gültigkeit dieser Beziehungen auch für 
spröde Körper vorauszusetzen. 

Man kann überlegen, daß unter dieser Voraussetzung J im Falle p,> p; stets negativ 
bleibt und deshalb bei genügend hohen Drucken der ‚plastische Verschiebungsbruch“ eintritt. 


3. Näherungsweise Berechnung einesersten Integrals 
a) Allgemeine Theorie. 


Beachtet man (14), so geht (17) in eine Differentialgleichung I. Ordnung über, wenn ent- 


weder o,, oder z, oder (1 — J) als unabhängige und |®| oder au als abhängige Ver- 


änderliche angesehen werden. Wir verwenden zund 


ee ie ae ir (22) 


und wollen nun versuchen, den Fehler bei einer Vernachlässigung von (1 + u) o, gegenüber |ß| 
abzuschätzen, da in vielen Fällen | |B| |> lo,| gilt. An die Stelle von (17) tritt dann ®): 


daR] = pet mit der Lösung: 


Y,.@= ar — Y„=const, das heißt nach (6): |S|=const . . . . (23). 
1—J 
4 d 
Die Differentialgleichung I. Ordnung läßt sich ohne Schwierigkeiten nach = na 


lösen”). (17) kann auch geschrieben werden®): 


(BIN _ arm. N, 


\14 ER ZER NE. 2 LAN, 
IB?“ £ 1 1+(ß —Az 
| =(4Wl): a an Eu er 
Ferner ergibt (14): 
ee. 
a De N 


2 er . 
so daß mit (22) und (24) o, jetzt als Funktion von n und z aufzufassen ist. (Für = >r et 
a © 
auch in allen späteren Gleichungen stets das obere und für = < 1 das untere Vorzeichen.) 


Aus (24) wird: 
2y2 n) 
er +An 


BE ee 
Für die folgenden Betrachtungen sei vorausgesetzt: 
k—-y2>0, und (BA) SS 0 nee een ne l27, 
so daß der Nenner von (26) stets positiv ausfällt. Aus (24) folgt: 
en} 
Ferner ist immer 
Al zu edchrauchh Agnes 07er Aare (29), 


dz 
unter den Voraussetzungen (27). 
Die Differentialgleichung (26) 


7(@)= 4Alz,n(); n(20) = No; 


6) Im folgenden mögen Näherungen durch Schlangenlinien bezeichnet werden. 
?) Die Ableitung nach z wird durch einen Punkt über der abgeleiteten Größe gekennzeichnet. 
8) Benutzt man nicht (6) sondern (1), so gelangt man zu: 


\ = z( % = )) lt 


En 
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werde durch eine Näherungsdifferentialgleichung: 

Y(z) = A(z, Y@)+ Da, Ye); Ya) = Yo: 0: (30). 
ersetzt und der Fehler zwischen ihren Lösungen abgeschätzt. Zu diesem Zwecke wurden pro- 
visorische Anfangsbedingungen gebildet. Für z, wählen wir den Mittelwert des Intervalles der 


Hooke schen Lösung 27(%): : 
a) (22(0)-+2u(1)). 


(26) wird sukzessive gelöst: 


„am +JACmıl)d; ma) Mm; 


Mit der nullten Näherung: n,(2) = Y(z) — (Y, — m) - 

Es läßt sich folgender Sat z°) beweisen: 

Es sei: |D(2, Y(z))| < öfür alle |? — 2,| < u. ea 

Ferner gelte: |n„(2) — no| < vfür allen >0O und |z — z,| < u. — Wenn in diesem Recht- 
eck |A(2,n„(2))|< A ist, kann gewählt werden v > Au. — 

Wird weiter vorausgesetzt, daß für alle 7®, n'® und zaus den Intervallen 

M®a— y|l<v bw M9—nl<v bzw k—2|<u 
die Lipschit z-Bedingung: 
14, n®) — Aa m) < U m® — n®); 

erfüllt ist, dann gilt die Fehlerabschätzung: 


|Y(2) — n(@)| <|Y — ml een (een 1) ac Lean! (31). 


Gelingt nun noch die Bestimmung von Y(z), «, Y, und n,, so kann die Fehlerabschätzung mit 
Hilfe dieses Satzes durchgeführt werden. 

Um bequemer rechnen zu können, werde die Voraussetzung der Lipschit z-Bedingung 
erweitert, indem nach dem Mittelwertsatz gesetzt wird 


M— max | | für alle 2 —z„|<u 
7 
und alle n im Bereich 
— v+-Min { Yo; <n <v+Max {Yı; nd. 
Im folgenden beabsichtigen wir drei verschiedene Näherungslösungen miteinander zu 
vergleichen, die durch verschiedene Indizes gekennzeichnet werden. 
b) Spezielle Fälle. 
Das ‚‚erste Integral“ (23) ist in der Tabelle 1 unter Ill eingetragen. 
Unter II wird die Näherungslösung Y;(z) behandelt, die sich bei Verwendung der Hooke- 
schen Differentialgleichung (19) ergibt. 
Y;(z) hat in den später betrachteten Fällen keine praktische Bedeutung und dient hier 
nur zum Vergleich. (Beachte Bild 3.) 
Die unter IT aufgeführte Näherungslösung Y,(z), welche aus dem Ansatz 
5, 
=) ==. 00 ART EEE 
i(2) Fr (32) 
hervorgeht, soll das erste Integral (23), in welches es für x — 0 übergeht, verbessern. 
Die Lösung von (32) lautet nach Y,;(z) aufgelöst: 


{ D ; (P —4 A) TE 4 
li fr 2) | . % (1—y 22)? a ill! BR = Re, Ar 
1 = : S 
Bee tu 
(33). 


ö, ist eine obere Schranke von 


A, (2, Y, (2)) = A(z, Y,(@)| — |D, (2, Y, (2))| <ö, 


Br Aa Yz@)l<d;. 
°) Kamke, Differentialgleichungen reeller Funktionen, Leipzig 1930, S. 94. 


und ö, von 
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Tabelle 1 (siehe auch Bild 3) 


Näherungsdifferentialgleichungen und deren Lösungen: 


ae ch ae N 
ar =) Be a a Van | 
er = TE EN 
zo. 
1 iR el, 
(explizite in (33)) 
| 87-0 
5 2y2 
Y= — 
2 (7 ir) Il 
ve a8, 
° 1—- y22+28,(ß—Az2) 
| III 
Ys=0 
Y,;=Yy 


Als Hauptfehlerglied erweist sich im folgenden jeweils der zweite Summand auf der 
rechten Seite von (31), und die Größe dieses Gliedes hängt im Wesentlichen von 6; ab. Man 
vergleiche dazu Bild 3. Dort sind die Y,;, nur wenig verschieden. Ebenso liegt n,, das auf Grund 
einer numerischen Methode gefunden wurde, zwischen den Y;,, so daß im folgenden nur die 
Fälle betrachtet werden sollen, in denen das erste Fehlerglied klein ist gegenüber dem zwei- 
ten Fehlerglied. Da, wie gleich gezeigt wird, für ein passend gewähltes x: ö6,<ö, ausfällt, hat 
unter dieser Voraussetzung Y,(z) im Vergleich zu Y;(z) einen kleineren Fehler gegenüber n(z), 
falls die zugehörigen Lipschit z-Konstanten sich nicht allzusehr unterscheiden. 


Mit z>0, o, <(0, haben in den beiden Fällen: 


1) ze < lıasmitseinemzei# 302 | 
: lade 


2) > 1 mit einem „>0 


A,(z,n) und A(z,n) (siehe (29)) gleiches Vorzeichen, wenn nur x nahe genug bei Null liegt. 
(Das ist ersichtlich, wenn der Nenner von A, (2, Y,) umgeschrieben wird: 


fl x [1 r) 0 MPer 42) 

at atz)+ re u=ym)) 
Für x = 0 folgt A, (z,n) = 0 und ö, = ö,. Übersteigt nun |x| nicht einen gewissen Wert, so 
fällt mit der Vorzeichenbedingung für x (34) 6, < ö, aus. 


emristisckes Vertahren. zur Bestimmung‘,von*% 
x sei so gewählt, daß Y;(2) und n'(z) sich nicht allzustark von einander unterscheiden. 
20 =) 
Es soll die Mittelbildung: JS (Yı(@) —n(2)) o(z) d& verschwinden. Wählt man speziell 
0) 


z2(i= 


f DT A 
die Gewichtsfunktion zu o(2) = ne benutzt (32) und (24), so liefert die jetzt elementar 


7 
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ausführbare Integration eine Bestimmungsgleichung für «: 
% 


(mr — [Bleu Zu Pa Eu 
Pi |Ple=0 — * Pi 
Verschwindet p, oder p;, so wird x — 0. Über deren weitere Verwendung siehe (41) und (46). 


4. Näherungslösung für die radiale Spannung 


a) Fehlerabschätzung bei der zweiten Integration. 

Um auch die Auswirkung des Fehlers zwischen Y;(z) und n(z) auf die endgültige Lösung 
o, (%) zu erkennen, wird dieser im folgenden in der Art berücksichtigt, daß die Rechnung auch 
Gültigkeit besitzt, wenn das nicht geschehen soll. Es sei gesetzt: 


9) - Ye) mM): - nenn: 6). 
Für die unbekannte Fehlerfunktion ;((p + 2) 07, o,) gilt die Abschätzung (31). 
Zur weiteren Behandlung wird der Fall x + 0 verwendet 20) 


ren G+Y(P+n) 0); 


Be 


2 ’ 0, O; 
v((p = %) Or; 07) >= 91((p + %) Or; 0,) + % == 
1-J J 
(Statt » könnte auch 9, geschrieben werden.) 


+ 3 -2 -7 0 1 2 3 Pr Is 


(el. 
Bild 2. (®) — — +1 und die Näherung 
r 


o, 1 g 
BER 2 als Funktion .von —*. 
„23 = 9% 


(37) ist — von y abgesehen — eine Differentialgleichung erster Ordnung für o,(x), die 
nach |) aufgelöst eine Differentialgleichung mit getrennten Variabeln liefert. Die Wurzel- 
ausdrücke für 2: und 2 erschweren diese Auflösung beträchtlich. Bei vorgegebenen Fehler- 
grenzen gelingt in gewissen Bereichen von o,/o, eine Näherungsdarstellung für 


eg. 


Bild 2 vergleicht die Wurzel mit 8 ZN \ In den Bereichen 


— 
| 2 


{02 [07 
—>2 bzw. — 
07 0, 


bleibt der Fehler kleiner als 15,5 %,. Man findet dann: 


Sl u, 


22 2 


BONN 


81-2, 


a. 0 ei \ i ER: 
) > ist eine Funktion von x. Da aber die mit x multiplizierte Differenz meistens vernachlässigbar 


ist, wird zin den Fehlerfunkti : a F Yeah 
Vo ine eriunktionen nicht als Argument erwähnt, zumal keine besondere Schwierigkeit durch deren 


. 
4 
+ 
& 


N N 
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Im folgenden teilen wir die Ausführung. 


= S—1d.b.:1) 1820: 


Hier kann das Glied mit Z unter der Voraussetzung (8) vernachlässigt werden. 2? ersetzen wir 
durch: 


2 —=|$l uote z((p+ 8) 07, 07) © 
(37) lautet nach |P] aufgelöst: 
BI=|Rol ++ Yl(p+ 2) or, o,) ; 
Bi: a hlore 
ER (2 be; 


} (1-—J 
Klp-+ 2) or 0) = IB ur ne 
Diese Lösung werde mit 7 = 0 zur Bestimmung von x nach (35) verwandt: 
pe\Erure 1-4 ym— rpulß+Y) 
= ————— I... .2.2202022. (A). 
pi 1-3ym—rm(P+Y) 


Nach der Voraussetzung (21) muß p,<p; sein, und damit fällt, was in (34) unter 1) voraus- 
gesetzt wurde, x <O aus, sofern, wie sich das in praktischen Fällen auch zeigt, 


3 
5 Y 
I1+u+xr>0 und —n<- 
; B+Y 
ist, 
Een 
07 0, 
Für 2? und z schreiben wir nach (39 
-- alı 57 ala + nab0)); 


&((p+ 8) or, 0). 
(37) lautet dann: 
(Rt BiRltr (RI 5) Hay (mie) 
D((p+ 2) 0, 0,); we, 
(+00, Dr NER pe | 


Sie stellt — von abgesehen — eine kubische Gleichung für Y— |B| dar. Wenn |o,| hinreichend 
klein ist gegenüber dem Absolutbetrage von |®|, läßt sich |'P| nach diesem Parameter ent- 


wickeln: x 
Bl! Haug, al) +) PERLE: . (43). 


Dies legt nahe, im ganzen Bereich Ro; (ähnlich wie im Falle ® < — 1) die Gl. (42) um- 
zuschreiben in: = Or 


IP] = |Pıl +0,0,+W((p+ 2) or, 0,) RE. or je (AAN: 


11) Die Zeichen < und > sollen andeuten, daß die Intervalle auf Grund vorgeschriebener Fehlergrenzen 
festgelegt sind. 
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— 


wobei versucht wird, jede Abhängigkeit vom linearen Gliede o, mitzubekommen. (42) ist 
erfüllt, wenn gesetzt wird: 


1 
8.) = ; 
Fe u au a 


BE ya TR) — at (7 a TR 
“nz ir. fa, 


1 AB Ro] 
Yp+ 2) 0.0) RO RN TE) 0 
+81 ANZ TEN) -VTR) o) 


m 


Hiermit sei im Falle: > 2; » graphisch berechnet aus: 
- > oo 3 — TE] 
* S. 2 23 Pe >: Ich 51 ie 
7 De sy as are ya Bl w) + (46). 
1 
Neil — zAlBol Y— |®ol 


zist hier noch eine Funktion von |®,| bzw. c,. Mit (21) erhält man positive Werte für x, wie in 
(34) unter 2) gefordert, wenn sich erweist: 


3 ir Re] Bol 
B—y— Zi Bl>0; x< a 
€ m Bol 


In den Bereichen (38) läßt sich die weitere Behandlung wieder gemeinsam vornehmen. 
Mit der Substitution: 


S=log(p+a); 
geht (40) bzw. (44) mit Hilfe von (14) über in: 
I d 
= Bl ( 2-04 2,0): 


Diese Differentialgleichung wird mit 


m 


= Pl - 2 — m) 2, 


verglichen und die Lösungen gegeneinander abgeschätzt. Zwei Anfangsbedingungen kann 
man dabei verwenden: 


für =0 folgt E= log p und (oe = (a (48a) 
oder 

für 1 folgt E=log(p +1) und (Jen = (Zen = — Mm - (48b). 
(Die nicht benutzte Randbedingung bestimmt die Integrationskonstante c, bzw. Y,..) 


Aus physikalischen Gründen wird vorausgesetzt, daß 0,(%) monoton zwischen —p; und 
— 9, verläuft, 


Weitere Voraussetzung: 
a ee 


A 
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Die unten beschrie- 


2) 


bene Lösungsmethode konvergiert, wenn in dem Rechteck: 


Die Grenzen des Intervalles ge sind abhängig von vr, 


o, zwischen — p; und — 9,; 


do, a 
de zwischen |® | +2 — a); FP und || +2 — %)Pa + P 0 
= 2 = 
1 Fr <P Bleibt. (für alleo, ee Par N). 
Daß = dann zwischen den N Grenzen liegt, kann man für die Bereiche > 2 
One 


(obere Vorzeichen) und >= — 1 (untere Vorzeichen) mit (21) und (49) leicht überlegen. 
Mit der een (68a) findet man o, wie folgt: 


do,(n don 
er | Pol— ee . N, nn); 
HA); (me) =— Di: 
Es gilt der Satz: ehrt. ni der Ungleichung (51) im Rechteck (50), so lautet die 
gesuchte Abschätzung: 
N r 
n— — cn PAR 
Ir a » 1 Sn 
Hätte man die Anfangsbedingung (48b) benutzt, so wäre man zu: 
= (e ai ar \ 
lo, — 2,|< en ne E Fer De 
gelangt. Berechnet man nun die Integrationskonstante aus: 
Bol 
ERERFRR TB -(? ren (54 
== Ze 
I Bol 19 p 
2 rw & . 


so werden die beiden auf Grund der Anfangsbedingungen verschiedenen Funktionen 2, in(52) 
und (53) identisch und erfüllen beide die Randbedingungen an den Stellen 2 =0 und z=1. 
lo, — 2,| ist für verschiedene # jeweils kleiner als das Minnimum der rechten Seiten von (52) 
und (53). Dieser hier angegebene Fehler hat seinen größten Wert an der Stelle <= Yp(p +1) 
— p <}. In der Nähe dieses Wertes von x aber verschwindet das Hauptfehlerglied von 9,, 
da dort 2 (&)= z,liegt. Das bedeutet: Der Fehler ist in der Gegend der Stelle & viel kleiner, 
als er hier angegeben wird. 

Die ‚endgültige Lösung‘ lautet: 


a (A Pre; 


2 — 2 — 0, »E 
(2 — ,) EIER ERERG, 
per ra re) 1; O<s<yp@tn— BL: - - 8) 
Z 5 ee 
0,(2) = Blade) +P(lp ER) 0a 100). 


Liegt «, zwischen 0 und 1, so hat lo,|, wie man sich in den Bereichen (38) mit (21) über- 
legen kann, seinen größten Wert am Innenrande des Rohres, wenn P genügend klein ist. 

Ebenso wird — von F(x) abgesehen — die Gültigkeit der oben gemachten Voraussetzung, 
o,(&)solle monoton verlaufen, durch (55) bestätigt. 


b) Numerische Methode zur Bestimmung der Integrationskonstanten c;. 


£ ö e ® x © > 102 . 
Die numerische Bestimmung sei zunächst für den Fall => 2 näher betrachtet. 
0, 
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Führt man nun folgende dimensionslosen Abkürzungen ein: er ; 3 
% a=2 — 0; Ben ER ee 


| I Bol 
h=Bm—-Pp); N=YM—Pp); A=Alpa — Pi)Y(Pa — Bi); an, 68); 
so heißt (54), (46) und (45): 
Bes) 
RE - (59); 
years Be 1 
Pi 2 
Pa Pa 
WS ee I Be 3, yo) 
| alt (nt zo), Ä a ıY® 
PR easy Dar 
2) I+u+% “ Di pi (60); 
p 5 1, 1 er 1 = 
+51 Pa (n+ 510) =, -( Yı 5 Yo) 
Pi pi 
1 = ee 
>(n0+5170°) +x(1+ß,0 — 1®—4,o®) 
ad : N Tv 


14. fo® 


Könnte aus diesen drei Gleichungen ® und x eliminiert werden, so bliebe eine Gleichung für a 
übrig, die von 6 dimensionslosen Parametern abhängt: 
y} Pa p + 1 3% Ru 
1 RE) TEE SZ enaae 5 W 
Dr p R 
Bei der Auswertung wird man zweckmäßig zunächst & aus (59) für verschiedene Werte 
von a berechnen. Zu jedem so gefundenen w gehört nach (60) ein x, und mit den bejden Werte- 
reihen ergeben sich aus (61) wiederum verschiedene a. Graphisch bzw. durch Interpolation 
bestimmt man dann ein alle Gleichungen erfüllendes « und nach der ersten Gl. aus (45) cı. 


(59), (60) und (61) gelten auch für Au — 1, wenn y, durch — y, und A, durch Null 
ersetzt wird. a 


c) Näherungslösung für « bei dünnwandigen Rohren. 
Die anschließenden Betrachtungen gelten in dem Bereich: 


u; ßı» Yı 


. 


R A 
ie 102: dh. we mp 210; wobei vorerst con Pe > - 


RE a 


sein soll 12). 
Zunächst wird untersucht, wie groß die dimensionslosen Parameter: 


= Ble+3): n=snlp+t2); »=zalp+z)Vlr+5) - + (63); 


ausfallen können. 
Die Misessche Fließbedingung: 


1/2 
(VER?) zruch — Vz lop| SEE NEE (64); 
[op] ist die Druckfestigkeit; verlangt bei verschwindender Wandstärke: 
1 


= bleibt auch bei nicht verschwindender Wandstärke erfüllt.) Das bedeutet nach (63) und 
Ye 


1 
R<Zpla]; n<5rlool; 1< > Alon|Yleo] RE (66). 


2) Es bleibt dann = 3 


er ee 
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Es läßt sich eine allgemeingültige Abschätzung finden, wenn man das Verhältnis. der 
Druckfestigkeit zur Zugfestigkeit kennt. Nach den Ausführungen zu (4) sei festgelegt: 


2 
ideal 3 
Hal Fun Be in 


a ist die Zugfestigkeit. 
Die Druckfestigkeit des Gußeisens ist nach den Be von Jung h lut hs) etwa 


viermal so groß wie seine Zugfestigkeit. Da es hier nur auf eine Abschätzung ankommt, sei 
für alle spröden Körper angenommen: 


Ion] s4 Ozug : 
Mit Hilfe der Ungleichungen (66) folgt dann: 


Pet Y2<7 
Eingangs wurde erwähnt, daß bisher stets B> C ausfiel, d. h. auch: 


2) Ba>»y.- 


Erfahrungsgemäß tritt beim linearen Druckversuch der Bruch bedeutend eher ein, Be Jden 
Wert In erreicht. (vgl. Bild 1). Zu In gehört nach Gl. EN 


: 
Yan] = a 


4 


Es sei vorausgesetzt: 
B 
lop! <z In ES A Re Re et re Lee ea (67). 


(Dann bleibt die erste Voraussetzung von (47) in dem hier betrachteten Bereich erfüllt, wie 
mit Hilfe von (72) erkennbar wird.) Die beiden letzten Gleichungen ergeben mit (65) 


m :V; (B—1)< 052 8, ; 


und es folgt als notwendige Bedingung: 


Al 2 
R<z: SON N EIER rt 
Führt man nun in (60) die Abkürzung: 
Pa 
Pi 
ein, logarithmiert beide Seiten, entwickelt die rechte in eine Potenzreihe von der Form: 
| a De Tee el 
— —— In (d+1)=1In ==) 
rare U 9 ur: ss +: 


N ee yo) | 
TERN LEE)EIMERD ER eerT 


und berücksichtigt nur das lineare Glied, so bleibt der Fehler kleiner als 0,8 %, denn mit den 
eben gefundenen Abschätzungen und (47) gilt: 


3f 1 = 3 A (0) 1 
= > 0)- Eh a NE 
<a ntzare)=5 br lernen SR 
wobei noch (72) zu beachten ist. Man findet so: 
1, 75 1 1 1 3 1 0) | # 3% 0) 
«(145 A10) L dr 79 a) n+3 Ale Zee Stute] 


A+u) [3 las at 
en amt za) [an Zara) 


13) Werkstoffhandbuch Stahl und Eisen, Düsseldorf 1937, L11, S. 5. 


= 
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Im Bereich (62) besteht die doppelte Ungleichung: 
1 1 1 1 ) 
Se Eee Ze 
3 >( Tanne 
so daß der Fehler bei einer Vernachlässigung des zweiten Summanden auf der linken ‚Seite 
gegenüber dem ersten kleiner als 10 % bleibt. (10 % beträgt der Fehler für d— x. In diesem 
Falle jedoch geht «— 0. Für d — 4 aber beträgt der Fehler höchstens 0,5 en! 
Mit der Abkürzung: 


De ._ u. ge De OB 
(p+3)1n (d+1) 


gilt die Näherungslösung: 
1138 1 0) 
o(p+4)(ntzAVa 
ER 1 Be 
hate role jan Far 


Wie man ohne weiteres sieht, zeigt sich die zweite Voraussetzung von (47) als erfüllt. 
Die rechte Seite der Gl. (59) kann in eine Potenzreihe nach 1/p entwickelt werden: 


1 were 1 
tz t gt en ie 


Vernachlässigt man alle Glieder mit 1/p und höheren Potenzen, so bleibt der Fehler kleiner als 
0,21%. Unter Benutzung von 


a 


-1 il 
a ee 
= 1% si25; 
P+Z er We) 
gilt: 
I 2 Be 
o=(p+3)@: o=1+at; Isusz 


==, 2 
Die Näherung Yo = 1+ = I — = t? ist mit einem Fehler kleiner als 1,05% und 


; 3 3 1 
Br = a ee 
Yo ie art get 
mit einem Fehler kleiner als 0,17 %, behaftet. 

Setzt man diese Näherungen in (61) ein und beachtet (63), so entsteht eine Gl. 4. Grades 
für a, die von den fünf dimensionslosen Parametern ße» Ya; As, D und t abhängt. Es scheint 
unzweckmäßig, diese exakt auszuwerten, da die brauchbare Wurzel schneller mit folgender 
sukzessiven Methode !*) zu finden ist: 

SER t | 3 ) t | t ) t?] 
a Ati —— — DI1+— 12 nn — — 
een Er 


; 
1+2Dß+2,[1 —{ ')a sp „or, 


ee 


r za, 
2(1+2Dß)+ 5m —9D)-3(1+2) 


U — 


1+2D&+4,[1 : ca 3D)+y9,(t—2D) 


| 14) 8 chu lz : Formelsammlung zur praktischen Mathematik, Sammlung Göschen, Bd. 1110, S.45 gibt 
eine Bedingung für die Konvergenz von (73) an, welche sicher erfüllt ist, wenn entweder 


1 
in U a5“ )g < 0,48 oder 
in 0<S%A< 0,7 I< 0,23 bleibt. 
Bei unverändertem Bereich des Radienverhältnisses bedeutet 1< 0,23 auch 2 >11,52. 
ı 
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Tabelle 2 gibt für u —= 0,3einen Überblick von der Größe des endgültigen Fehlers bei der 
Berechnung von a. Der Fehler zwischen den zugehörigen Werten bleibt in der Tabelle kleiner 
als 0,22 %. 


Tabelle 2 Einige Werte der Potenz a nach der Näherungslösung (73) (obere 


Ziffer) und den exakten Bestimmungsgleichungen (59) (60) und (6l) (untere 
Ziffer) für u=03. 
1 . = 
= m 0,05 0,1 
2 p 
Fr 2 0,2 0,4 Fr Y 0,2 | 0,4 
B,=0,6 1,351 301 0,843 494 E06 1,374 628 0,895 347 
EYE 1,351 158 0,843 263 =, 1,374 402 0,894 720 
0,1 0,1 
B, = 0,8 1,376 035 0,887 442 — 08 1,415 148 0,967 268 
2 ©°|| 1,375870 | 0,887441 P2=0,8|| [’414851 0,966 447 
1,5 en 
086 1,270 595 0,795 115 er 1,283 646 0,827 877 
z : 1,270 503 0,794 970 B=V0, 1,283 792 0,827 594 
0,2 EEE STEIEREEE AZ ger 0,2 
B. = 0,8 1,296 604 0,838 610 08 1,327 665 0,901 243 
a : 1,296 502 0,833 368 Da 1,327 735 0,900 743 
(t = 0,097.561.0; D = 0,124.141.4) (t = 0,190.476.2; D = 0,242.371.3) 
hl 0,2 0,4 = Ya 0,2 0,4 
B,= 0,6 1,372 682 0,888 755 — 06 1,403 656 0,960 055 
2 2 1,372 880 0,8838 668 R=, 1,405 083 0,960 289 
0,1 0,1 
B. = 0,8 1,410 622 0,956 337 RB, 08 1,460 607 1,059 977 
2 ; 1,410 800 0,956 323 2 ©” ||. 1,461 994 1,060 200 
3 
B.= 0,6 1,283 581 0,824 397 B = 0,6 1,298 236 0,868 349 
2 ; 1,284 030 0,824 471 2 ? 1,301 069 0,869 311 
0,2 0,2 
Bs = 0,8 1,324 594 0,892 927 Br = 0,8 1,362 924 0,975 141 
= i 1,325 072 0,892 985 E 2 1,365 661 0,975 300 
(t = 0,170.731.7; D = 0,220.432.9) (t = 1/3; D = 0,430.369.0) 
Wird in (73) a als Funktion von 1/p betrachtet, so liegt der kleinste Wert von a stets an 
der Stelle ——=0, d.h. auch = 0 und D=(, wie es sich numerisch zeigt und bereits aus 


Tabelle 2 zu entnehmen ist. Es gilt dann: 


(za) 


Die rechte Seite dieser Gleichung aber ist größer als 0,2. Daraus läßt sich schließen, daß in 
allen Fällen 2 wahrscheinlich «> 0,2 bleibt, womit die Voraussetzung (49) erfüllt ist. 
[07 ae] 


5 
Die Grenzformel für kleine Wandstärken (74) läßt sich als „Faustregel“ bei dünnwandigen 
Rohren verwenden, denn Tabelle 2 lehrt, daß a nur geringfügig von 1/p abhängt. 


Analog findet man in dem Bereich ®): 


a! 
nen... 


R 
Is <s1ll; = 
Pa 


=R, 


(6) 


Oo 
15) Es gilt dann stets a —1. 
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die Näherungslösung: 


2(1+2Dß)+33[143) 
A eeriernt 2 Rei vi): 
Vene; 02%>—0,15; ir 0=2D> 5 


und die Grenzformel für kleine Wandstärken: 


[e=2+3%] 9 
d) Zur Bestimmung von ng. 


Bei verschwindender Wandstärke streben die Fehlerfunktiorfen in der Lösung (55), (56) 
gegen Null. Je größer die Radienverhältnisse ausfallen, umso größer wird der Fehler, und man 
kann dann versuchen, den quadratischen Ansatz von (43) zu lösen, was auch elementar durch- 
führbar ist, wobei sich aber der Fehler nicht wesentlich verringert. Mit Hilfe dieser Lösung 
gelingt es jedoch, die Differenz Y, — n,, die man zur Fehlerabschätzung (31) noch kennen muß, 
ungefähr zu berechnen und notfalls auch abzuschätzen. 

Im folgenden sind lediglich die Formeln wiedergegeben, mit welchen eine ungefähre 


Bestimmung von Y, — n, möglich ist. Sie gelten für a 2. (Im Falle =< —1 mußy 
durch — y, A durch Null und y—|®,| durch Y|R = ersetzt FE ) 


2m "ey, m+W et 


en Bl : 
play I mW 
Pol | 
ey 
W=(2—%) (Pat Pi) + (Pa — Pi) ES a: 
we 


a=y2 — 9% -—40,; 


Bel RP: eo 
%,=(2PQS— PT)(%) — RP%; 


RT TR 
‚sel L q\ 
= ar Re at) 

1 
m en 3 4 ni j 
"373 v8, (t 204 2) 5 


P=1—(B— y)|Bol+2 [Bol Y— Bol; 
= la TR) + TR tr. 


i und q kennzeichnen die Annäherung von D(z, Y (z)) durch eine Gerade: 


Da, Y(@)=—i@-2)+g. 


Zu dieser Annäherung beachte man Bild 6a und 6b. 


Mit Hilfe dieser Bestimmungsmethode läßt sich wenigstens nachprüf 
t Hi en, ob das erst 
Fehlerglied in (31) gegenüber dem zweiten vernachlässigbar ist. i Eee 
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e) Durchrechnung eines Beispiels. 
Nach diesem Schema wurde ein Beispiel behandelt. 


22,5. 10° kK8r cn, u—=0,3;; E =.10C kg cn 2 en (77). 
R,\? 
a =1,3; 2a = 000 kgcmz2; pP: == 100. kg? 
(7 


Der Zustand des Rohres liegt in der Nähe des ‚plastischen Verschiebungsbruches‘“. 


Folgende Ergebnisse seien festgehalten: Sowohl o, als auch Y(P’)? verteilen sich gleich- 
mäßiger über den Querschnitt des Rohres als die entsprechende Hooke sche Lösung. 


7@ S 
| N 66 
| 52 
| 
I 
-7720 | 
j 
| 
I 
| 
| 
| 
j N 
-7740 | ee 
| al 
| | 
Kress | 
L | 
| 
-1760 | 
| 
| 
| 
| 
N | 
5 | 
-7780 N l 
N » | 
N, 
\ Ä a 


Bild 3. Die ersten Integrale mit ihren Fehlergrenzen für das Beispiel (77). 


Die in Tabelle 1 eingetragenen drei Näherungslösungen wurden an diesem Beispiel 
erprobt, um zu sehen, in-wie-weit sich die Fehler der Y;(z) gegenüber »(z) nach (31) unter- 
scheiden. Sie sind in Bild 3 schraffiert. (Dabei wurde der Übersicht halber jeweils das Fehler- 
glied mit dem Faktor |Y;, — no| vernachlässigt. Sein Wert beträgt an den Grenzen der ver- 
schiedenen z-Intervalle nur rd. 1/20 des gesamten Fehlers.) 

Y,(z) und Y,(z) bedeuten numerische Löusngen des nächsten Abschnittes. 


f) Brauchbarkeit der angegebenen Lösung und allgemeine Ergebnisse. 
Mit Hilfe der Näherungslösung aus (55), (56): 


— & \ 
(8) = — sl 2 ); S=N-+ 9; 


| 

p+% [ N ‚| / 

0,(2)= En Ss SE 18) 3 S = (et N (? 9 ) | Te (78); 
p 


läßt sich nun angenähert bestimmen, welche Bereiche des Druckverhältnisses durch (38) aus- 
geschlossen werden, wenn man a als Parameter betrachtet. Es ergeben sich folgende Bedin- 
gungen: 

1. Neben 


(2-@e-o(?#2\)>o 


\ 


-I 
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muß Ki 
Fr) 
PiS p 


"le Re 
I 


sein, wenn 2 < — 1 ausfallen soll. 
0, 


2. Neben a 

(er 2) )> 
p 
muß 
Pa a N 
z z 
Pi “ -) 
a 
T p 


; 0 
sein, wenn => 2 ausfallen soll. 
0, 


Die Voraussetzung (21) ist in der Näherung (78) für alle a erfüllt. 

Die Voraussetzungen (27) sollen jetzt auf ihre Brauchbarkeit untersucht werden. Im 
Bereich des ‚‚spröden Bruches‘“ mit |B]> O gilt mit J < 0,72 sofort 

1—-yz2>0. 

. Beim ‚plastischen Verschiebungsbruch‘“ verlangt (64): 2? <|o)|. yz? kann also nach (66) 
höchstens die obere Schranke von 2 y, erreichen. Diese ist in (68) aber 1,2. Das nicht-lineare 
Verformungsgesetz (6), das das spröde Verhalten eines Körpers wiedergibt, besitzt aber nicht bis 
zum Eintritt des ‚‚plastischen Verschiebungsbruches‘ Gültigkeit, da bereits vorher mit pla- 
stischen Vorgängen gerechnet werden muß. Wird die Voraussetzung gemacht, daß der Gültig- 
keitsbereich von (6) durch 


5 


beschränkt ist, so bleibt sicher y, < 0,5 und es fällt 1— yz2?> O aus. (Das ist auch eine zu- 
sätzliche Begründung von (67).) (B—A2)> Oist auf Grund von (64) und (67) stets erfüllt. 

Damit sind alle während dieser Arbeit gemachten Voraussetzungen auf ihre Brauchbar- 
keit hin untersucht, und es konnte wenigstens in speziellen Bereichen ihre Gültigkeit nach- 
gewiesen werden. Als eigentliche Voraussetzungen bleiben lediglich 0, =0, o, SO und 
(1— J)> 0,28 bestehen. 

Es läßt sich leicht noch näheres aussagen über die Größenordnung des zur Lösung (55), 
(56) gehörenden Fehlers, worauf jedoch hier verzichtet werde. 


Unter welchen Bedingungen verläuft nun /(P’)® nach dem Ansatz (78) monoton über den 


Querschnitt mit einem Maximum am Innenrande, d.h., wann ist nn. <0 ? Mit (14) 
: d 
folgt im Falle 9,< p; ist stets Fr <0; und im Falle 2,> p; nur dann <0; wenn 
gleichzeitig die Beziehungen gelten: 
el, 
Pa ı 
a>—; al > — 
2 aa m; f a 
( ) Di le, Ba 
Tee “ ” 
a2a—]1) 
Mit den bisher gemachten Voraussetzungen läßt sich nicht allgemein zeigen, daß stets a <0 
x 
ausfällt, m nach (73) a bis 0,2 abfallen kann. (Verlangt man (79), so bleibt allerdings «> 0,5%) 
Ist ER <.0, so folgt auf Grund der Beziehungen (28) und (29): 
dn d|6© d(1— J) 
7) i d.h. nach (6) auch: ac und wel 


ee Bee 
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Die Maxima von |©| und (1 -— J)liegen 
am Außenrande des Rohres. Im Falle 
Pa <Pp; hat also die Anstrengung J stets 
an der Steller = R; ihren größten Wert, 
so daß bei genügend hohem Drucke der 
„spröde Bruch“ am Innenrande aus- 
gelöst wird. 

Abb. 4 zeigt den Spannungsver- 
lauf nach (78) im Falle 

Ra _ 9. Pa 
R, 2; 2 
für verschiedene Werte von a. 

Man kann auch a in roher Weise 
verhältnismäßig einfach ermitteln, wenn 
man (78) als Ansatz in die Differential- 
gleichung (17) einsetzt. 


5. Numerischer Lösungs- 
weg 

In diesem Abschnitt wird eine 
numerische Methode entwickelt, der 
eine sukzessive Approximation mit 0,4 
und o,z als nullter Näherung überlagert 
ist. Es sei hier die nach |'P]’ aufgelöste 
Beziehung: 


Br _ 
P+ © 
3(B| — 2 0,)?- Ei, o,, 0,) 


 (@+2)%(1+0]B]-3 6,): E(a. 0, 6,)) 


== J(%, 0,05); 
E(«, 0, ; 07) >= N 


( B| er (1 as m) o,) (7 E a Bild 4. Die Spannungsveıteilung im EprOgen Bob: Er Ansatz (78) 
22 £ für Potenzen verschiedener Größe, m RR - =2 
22 ye— (+) A2))' ee 


benutzt, die sich aus (17) ergibt. Sie enthält auf der rechten Seite o; nicht mehr. Die Sub- 
stitution: 


9,(%) — 3(2) +&(8) : 


homogenisiert die Randbedingungen zu: 


CO) TE 
& (x) erfüllt die Differentialgleichung: 
N Se BER 
+ p+s = 44,0; 65 


9(2,%,C) = I(8, 9,u(8&) +8, ul) +). 


Zur sukzessiven Lösungsmethode !%): 


1 
ea “% =0(0; Enti = 7% 2) ° 9(& u (1); Cr (21) da; . 
: N 
En+ı = es ; 91, &n(&ı)> &n (2) da, ; 
{0 
Die ezakte Integration der ersten Näherung erwies sich als kompliziert. Die Integrale können dabei auf 
elliptische Integrale 3. Gattung geführt werden. 


IE 
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gehört die Greensche Funktion (Einflußfunktion): 


E ern. 2 ware 


G(8, 2) = 2 für TE 
ee -eENerm: »anet 


Die einzelnen Iterationsschritte lassen sich leicht numerisch auswerten !”). Dabei zeigt sich in 
dem Beispiel (77), daß die erste Näherung durch die zweite nur geringfügig verändert wird, 
Eine Fehlerabschätzung dieser sukzessiven Methode läßt sich mit Hilfe eines Konvergenz- 
beweises von Picard!) führen. 


Bildet man hiermit wieder a 7 als Funktion von z, so ergeben sich drei Kurven, 


die ebenfalls in Bild 3 eingetragen sind: Zur nullten Näherung gehört Y, (z) und zur ersten und 
zweiten Y, (z) bzw. Y,; (2). 


12 Kapitel 
Näherungslösung einer Differentialgleichung für die radiale Verschiebung 


1. Die Abhängigkeit des Spannungstensors vom Verzerrungs- 
tensorund dessen Invarianten 


Aus dem nicht-linearen Verformungsgesetz 19) I (6) folgt: 


Da £ 2d 
B- (84, u DE ee 
1 — J kann geschrieben werden: 
BE ae 7 
le J) Dee 
LE“: B 


ISIY Sy — T° 


“r. 


1—2afita 
Das ist eine kubische Gleichung für Y(1 — J): 


| VE NEMIYA— NN 0 ge 
M=LE@%: N-1H2lelta@; 
SUTBELN 10 ONE E a 
er ae FR be 2 Iru | 
die im Falle Z—= 0 durch: a 
1 
a 
gelöst wird. Wenn man in der aus (2) folgenden Beziehung: 


= Ä 
N Myi 


den bei M stehenden Faktor Y(1 — J) als Konstante Y(1 — J), betrachtet und geeignet wählt, 
erhält man eine gute Näherungslösung: (1 — J). 


1 


In Bild5 ist 
1 PEN 
N1—J) = mit N1—J)= L 
De lee 
m} -ZÜU=7), 
verglichen, wobei die rechten Seiten vorher noch mit der aus (2) gewonnenen Gl.: 


I k . er TEN 4 ——— 
vll J=775; mit = MI) LIBRI) 


17) Vergl. v. Sanden: Praxis der Differentialgleichungen, Berlin 1945, S. 89. 
18) Trait6 d’Analyse. Paris 1905, Tome III, Chapitre V, S. 90. 
19) I( ) bedeutet die jeweilige Gleichung des 1. Kapitels. 
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umgeformt werden in: 
Nıenerıe NAZN=— — —. 
1 3[1- YO) 
Der Parameter a I) wurde zu 0,8 bzw. 12 gewählt. Es zeigt sich, daß für diese 


1—J 
extremen Werte (1 — J) im Intervall 0 > k> — 1,1 noch eine ausreichende Näherung ist. 


Bild5. Zur Veranschaulichung der Näherungen für N1—J)=1+k 


Im Hohlzylinderproblem sei (1 — J), durch die Hooke sche Lösung bestimmt: 
1 v 
Ü-NM=zi-In0)+ Jul): 
Man erhält dann für das Beispiel I (77) 
A—J) 1-Ih 
la a 
1-70) 022 und u) 
und der Fehler zwischen 1 — J und er) beträgt für k = — 2 erst 4,5%. 
Im folgenden wird die näherungsweise Umkehrung 
1 
1+J, 


—=.0,08% 


1—-J= 


benutzt. Dabei ist 


nr 
mn 
> 
— 


J, = B|&| + NS) — LISIY(SYAa — I) 


eine Invariante des Verzerrungszustandes. 


2 Die Differentialgleiehnng für die radiale Verschiebnng und 
deren Lösungsmethode 
In Zylinderkoordinaten ergibt (1) mit (4) unter Verwendung von I (12) und I (13) neben 
I (16) noch die beiden Beziehungen: 


De E 


E , Y Y ’ 1 
2 ee rer Greene er A 
Ferner gilt ?°): | 
— 2u 
Cpe ey) nn 0 


20) en darf nicht mit dem & des I. Kapitels verwechselt werden, 
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ee 2 
(= 3° : 


SEIT ea 
EEE LEE at 1) (7). 


oder 
4,7 1+u\? os BER, Yy ; Bine, Y 
o-\/e( 9 Bo : ee ı 2=y p+ 


(5) und die Gleichgewichtsbedingung I (9) liefern die gesuchte Differentialgleichung: 


: , Galler di) af’): 1+J), =1+bA+co!— 14. . (8); 
BE A wen: I 


3z = en 1— u 3 


ee. 
Die Lösung dieser Differentialgleichung muß nach (5) mit (4) den Randbedingungen genügen: 


Pi =} E 
os: 


oe, a a 
(8) ist dann eine Differentialgleichung I. Ordnung, wenn man z. B. A als unabhängige und ® 
bzw. Q als abhängige Veränderliche einführt. 


Hier soll jedoch an die Ergebnisse des I. Kapitels angeknüpft werden. Das liegt nahe, da 
die Differentialgleichungen (8) und I (24) identisch sind. Die dortigen Näherungsdifferential- 
gleichungen liefern: 


ee ee Zu a 
ee ut u )=0.... 1); 


h (E / en ’ x ‚ Ya Reg 
: () +) werk + ) nn 


In diesem Kapitel braucht man sich nicht mit deren Lösungen zu begnügen, sondern kann Y, 
als nullte Näherung bei einer sukzessiven Approximation der Differentialgleichung (8) ver- 
wenden. Dabei muß jede Näherung den Randbedingungen (10) genügen. Der prozentuale 
Fehler zwischen der nullten Näherung A, —=const und A liegt in der Größenordnung des 
Fehlers aus I (31) zwischen Y, und n. Bedeutet y, (2) die n-te Näherung und 


An, = ym+ ns 3 h2(&) = h(®, Yn > Yn, Yn) 3 
so laute: 
An = n-1(8) ; Mel... 2 u N (13); 


Die Lösung der nullten Näherung ist: 


D, 


ee Te seit); 
die der n-ten Näherung: 
Y% = Hn-ı(2) + D 1: Den+z . 
In n—ı\% an+ı (P + %) -| er ; 
ei / D n 
: ee 2). 


IH = 1 : 
p+2 =; | In-ı (8) de + | (?+ %)° In-ı(8) da 
2 2 . 
Bi_,(a)) Do 


Hn-ı(%) | 


Er Zr" 
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3.Bestimmung der In- 
tegrationskonstanten 


Liegen die Funktionen H,-ı (=) 
und Ha-ı (8) fest, so müssen noch 
die Integrationskonstanten Dan+ı 
und Dg„+2 bestimmt werden. Mit 
der aus (5) und (7) folgenden Bezie- 
hung: 


lu 
Br ae 2 = 
Be 


2, 


kann 1 — Jals Funktion von o, und 
x angesehen werden und ist an den 
Stellen &=0 und 2 =1bis auf die 
Integrationskonstanten bekannt. Die 
Randbedingungen (10)ergeben damit 
zwei algebraische Bestimmungsglei- 
chungen 6. Grades in Dan+1ı und 
Dan+2. Verhältnismäßig wenig Re- 
chenaufwand erfordert eine sukzes- 
sive Bestimmung der Integrations- 
konstanten, bei welcher jeweils 
1— J (0) und 1 — J (1) mit der vor- 
herigen Näherung berechnet werden. 
Die nullte Näherung soll dabei die 
Hookesche Lösung sein. In dem 
Beispiel des I. Kapitels ist 1 — Jz 
nur bis zu einigen Prozenten von 
1—- J eines spröden Körpers ent- 
fernt, weshalb die Näherungsschritte 
dort schnell konvergieren. Wir wol- 
len uns hier auf ein solches Verhalten 
beschränken. Eine Konvergenzbe- 
trachtung gelingt im speziellen Fall 
nach »Schulz?%): 


-A,(Z.% (Z)) 


zZ 
Bild 6a. Näherungsdarstellung für / D, ( Y, (z)) dr 
20 


2 
Bild 6b. Näherungsdarstellung für / D, (7, F, (0) dr 
2o 


für das Beispiel (77) 


Die dann für jedes feste n entstehende Folge sei ausführlich hingeschrieben, dabei steht 


(2) für ——— deren 


rd: 
| ) 
om(PrIY _« 
SEN a 
2n+1l — E an (p+1)? E ee 3 
p » 
v Pa E ’ (1) a 
ö ee Rs 17 LEE En ) (16) 
ml 2 Pi E ' | Hn-ı(0) N 
per — A — 70m | +u r ) 
1+3(0 (a)? 
— ZN 9% 9,(% ) 4 / 
(1— J(x)) Pa oe) +7) 7 A Brad) 
9, o,(&)1/1+ 3(0(z))?. 
1+ 0%(8) Yr (14 09 (a)? 3 ) 


21) Vgl. Anmerkung #), dort S. 62, 
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De : | Hn-ı(%) 
AIR 
ee 1—u (p+@)% n-1(2) + Fr 
a 0) H Hn-ı(@ 
2 Dont Hn-ı(a en eig 
Die Folge beginnt, » — 0, mit: (16) 
p+1Y | 
TE a ) 
aa Me 
El 
Q m (2? 
Dabei gilt: 0,(0) = — 9; und o, (1) = — M- N 


Bei der Bestimmung der Integrationskonstanten der nullten Näherung (14) verschwinden 
die Funktionen HZ, (x) =0 und Hl, (e) =. 


4. Näherungslösungen der Differentialgleichungen für die 
radiale Verschiebung 


Neben (14) kann auch die Lösung von (12): 
1 1+x 
vn) + er ee 
als Näherung verwandt werden. Dabei ist nach I (34) und I (21): 
1j’0, =. 1° turen, m 
2), Hr lsslürgn en 
Für o, galt die Näherungslösung: 


. (17); 


,=—8,+8 (2) mit 52. 


Im Hooke schen Fallist a, =1und a=2. Indem Gültigkeitsbereich der Näherung I (78) 
fällt |l — a,| bedeutend kleiner als 2— aaus. D.h.: Die Art der Abhängigkeit vom Radius r 
für die radiale Verschiebung s unterscheidet sich weniger von der Hookeschen Art, als es 
für die radiale Spannung o, zutrifft. Trotzdem aber liegt y von Yz prozentual weiter entfernt, 
als o, von 0,7, denn für o, und o,z gelten die gleichen, während zu yandere Randbedingungen 
als zu Yz gehören, was sich bei den Werten der Integrationskonstanten bemerkbar macht. 

Die Auswertung des 1. Schritts von (15) gelingt mit elementaren Funktionen, wenn man 
für & eine Näherungsdarstellung verwendet. Bei einem Fehler kleiner als 15,5 % findet man: 


Y 7 ——— x 
17 we Beste 1A: 
| ea u+uR; Dr 


w* Y . 
Ten me, 093. 
| D 4 2 Y u + u 0% nm ‚> 
Auf die Wiedergabe der Formeln für y, möge hier verzichtet werden. 
Die Werte von |S| am Außen- und Innenrande unterscheiden sich für das Beispiel I (77) 


bei der Lösung %, um 1,1%, bei %, um 0,8%. Bei verschwindend kleiner Wandstärke ns 1) 


bleibt exakt: |&|—=const. Und wie das Beispiel zeigt, kann man in guter Näherung auch bei 
nicht verschwindender Wandstärke diese Lösung dd. RE) verwenden. 


Eingegangen am 31. Januar 1953. 


F 
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Über erzwungene nicht=lineare Schwingungen hoher Erreger: 
frequenz und ihre Stabilität 
Von K. L. Stellmacher in Göttingen 


Um Binsicht zu gewinnen in das Verhalten von nicht-linearen Schwingern bei nicht-linearer periodischer 
Erregung hoher Frequenz, werden periodische Lösungen gewisser Differentialgleichungen zweiter Ordnung 
konstruiert, und die Stabilität dieser periodischen Lösungen weiter untersucht. Es ergeben sich gewisse Ver- 
en der periodischen Lösungen, die hier im allgemeinen gleichzeitig Stabilitätsgrenzen dar- 
stellen. 


Aim of the paper is to study the behaviour of non-linear oscillators in the case of non-linear excilalion o 

‚ high frequency For this purpose, periodical solutions of certain differential equations ofthe second order are 

constructed, and the stability of these solutions is investigated. T’he periodical solutions have certain points of 
bifurcation, that represent, in this case, generally limits of stability. 


Pour &tudier la reaction d’oscillateurs non lineaires dans le cas d’excitation periodique non lineaire de 
haute frequence, des r&solutions periodiques de certaines Equations differentielles de second ordre sont con- 
struites, et la stabilite de ces r&solutions periodiques est ewaminee. De certains points de bifurcation des 
resolulions periodiques sont Irouv£s, lesquels, dans ce cas, reprösent des limites de stabililE en general. 


B cTaTbe HsyyasTcH HoOBeNeHHe HEJIHHEÜHBIX OCHMINATOPOB IPH HeJIUHEÄHOM TEePHOAHYECKOM 
BOSÖY:KACHHH BBICOKOÄH YACTOTEI. Ü 3TOÄ MenbIO KOHCTPYHpymTcA mepmoNnyeckue permeHnusn 
HeKOTOPEIX AußhpepeHnuHanbHBIx YpaBHeHHÄa BTOPOTO HOPANKA, a 3ATeM Hecienyerca YcToä- 
YHBOCTb 9THX IIEPHONHYecKuX peMeHü. B pesyıIbBTaTe IONYYAMTCA H3BeCTHLE TOUKEH BSTBJICHHA 
TePHONUYecKuX pemieHHnf, ONHOBPeMEHHO IIPeACcTaBıAmımHe CODoH, B NAHHOM CIIY4Yae, IIPeAeJIst 
YCTOAYHBOCTH. 


Nicht-lineare Schwinger, auf die periodische Kräfte wirken, deren Frequenz einiger- 
maßen groß ist gegenüber den Eigenfrequenzen des Schwingers, zeigen unter Umständen 
starke Auslenkungen der Mittellage ihrer Bewegung aus ihrer Gleichgewichtslage. Ein Schul- 
beispiel dafür ist das Verhalten eines schweren Pendels, dessen Aufhängepunkt periodisch 
bewegt wird. Technisch wichtige Anwendungen, die hierher gehören, sind u. a.: Auslenkungen 
von Gerätezeigern infolge von Erschütterungen sowie Schlingerfehler von Magnet- und 
Kreiselkompaß. Alle diese Systeme sind dadurch gekennzeichnet, daß die periodische Er- 
regung nicht-linear mit dem Schwinger gekoppelt ist. 

Interesse verdient die Behandlung der Stabilitätsverhältnisse. Bekanntlich kann aus 
einer stabilen Gleichgewichtslage durch Erschütterung eine instabile werden und umgekehrt. 
Für die meisten der technisch interessierenden hierhergehörenden Erscheinungen sind bisher 
nur mehr oder weniger grobe Näherungen für die Gleichgewichtslagen veröffentlicht. Schwierig- 
keiten und Unzuträglichkeiten sind jedoch bei der Diskussion der Stabilitätsverhältnisse auf- 
getreten und zwar dadurch, daß die hier interessierenden Schwinger in zu starker Lineari- 
sierung von vornherein als ‚‚rheolineare Systeme‘ behandelt wurden. Tatsächlich sind jedoch 
viele der in diesem Rahmen interessierenden Schwingungsvorgänge so stark nicht-linear, daß 
es mir notwendig erscheint, auf die exakten Bewegungsgleichungen zurückzugreifen. 

Ich gehe im Folgenden so vor, daß ich zunächst eine möglichst allgemeine zeitlich peri- 
odische Bewegungsgleichung mit einem Freiheitsgrad zugrundelege und von dieser periodische 
Lösungen konstruiere. Ich definiere dann als ‚‚,Gleichgewichtslage‘ in naheliegender Weise das 
zeitunabhängige Glied der sich ergebenden Fourierentwicklung der gewonnenen peri- 
odischen Lösung. 

Die Stabilität dieser Gleichgewichtslage wird dann definiert als (infinitesimale) Stabilität 
der gewonnenen periodischen Lösung, indem untersucht wird, ob eine zur gewonnenen Lösung 
benachbarte Lösung für alle Zeiten benachbart bleibt oder nicht. So wird man in bekannter 
Weise auf die Diskussion einer gewissen linearen Differentialgleichung mit periodischen Koeffi- 
zienten (also vom Hillschen Typ) geführt, die sogenannte Jakobische Derivierte der 
vorgelegten Bewegungsgleichung. Deren Diskussion erfolgt hier in einer von der üblichen Art 
abweichenden Behandlungsweise, die der mathematischen Technik angepaßt ist, mit der die 
zunächst zu gewinnende periodische Lösung konstruiert ist, nämlich durch Entwicklung nach 
einem Störungsparameter. 

Es werden so gewisse mathematische Sätze allgemeinerer Art über periodische Lösungen 
von Differentialgleichungen gewonnen, die dann im letzten Abschnitt als Beispiel auf den Fall 
des schweren Pendels mit einem Aufhängepunkt, der auf verschiedenartige Weise periodisch 
bewegt wird, angewandt werden. Damit werden dann auch die umstrittenen Stabilitätsfragen 
für verschiedene Anordnungen eines solchen erschütterten Pendels geklärt. 

Die Untersuchung ist jedoch allgemeiner angelegt, als für die speziellen Beispiele des 
Abschnitt 5 notwendig wäre. Dadurch tritt das Grundsätzliche klarer hervor, und es sind 
eine Reihe ähnlicher Anordnungen mit erfaßt. 
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1. Konstruktion einer periodischen Lösung 
Wir gehen aus von einer Differenzialgleichung der Form 
d2y 
| Fr?) 
F* sei regulär analytisch in allen drei Veränderlichen, außerdem in 7 periodisch von der Kreis- 
frequenz zwei. Der „‚Störungsparameter‘“ a sei von z und der gesuchten Funktion @ unab- 
hängig. Darüber hinaus gestatte F* eine Entwicklung 
oo —+k : 
Bar 3. 2 area 
k=1l=—k 
(Den Index k nennen wir die Ordnung, I die Vielfachheit der Frequenz.) Wir setzen also voraus, 
daß die Vielfachheit der Frequenz bei keinem Gliede größer ist als die Ordnung k. Ein Beispiel 
ist etwaa F*—=G@(acos2r,«,_). In vielen Fällen der Anwendung wird es sinnvoll sein, eine 
solche Abhängigkeit von einem Parameter künstlich herzustellen, wobei es besonders auf die 
Anpassung an die physikalischen Gegebenheiten ankommt. (Beispiele dafür im Abschnitt 5). 
Im Folgenden wollen wir von der willkürlichen Vorstellung ausgehen, daß der Parameter a der 
Erregerfrequenz proportional sei. 

Die Tatsache, daß wir uns allgemein für Bewegungen hoher Frequenz interessieren wollen, 
kommt dann in der Gestalt der Differentialgleichung darin zum Ausdruck, daß kein von & 
freies Glied auftritt. 

Ziel ist zunächst, eine periodische Lösung zu konstruieren durch ‚‚Störungsentwicklung‘, 
d.h. durch Potenzreihenentwicklung nach dem Parameter «in der Umgebung von a =. 

Zu dem Zweck führen wir mit Siegel!) die neuen unabhängigen Veränderlichen ein: 

er. gl 2, ee rolg en: 0 A 
und die Bezeichnung ee ee .... 
Wir fassen dann alle zu bestimmenden Größen als Funktionen der beiden unabhängigen Ver- 
änderlichen & und n auf. Wir setzen fest, daß bei der Entwicklung einer Funktion nach 


Potenzen von &und n = 
® == % 1:05) 5 Apı &kı ELIE e  E Pe  r (1.3) 
k=1 |I|<k 


grundsätzlich alle Koeffizienten a;; gleich Null zu setzen sind, für die ||| > kist, auch ohne daß 
dies ausdrücklich dabeigeschrieben wird. Da diese Vorschrift bei den elementaren Rechen- 
operationen und bei Differenzieren nach 7 erhalten bleibt, und da wir schließlich keine anders- 
artigen Reihenentwicklungen als Ausgangsfunktionen zulassen, können dabei keine Wider- 
sprüche auftreten. Wir schreiben noch einmal in dieser Sie gelschen Schreibweise unserer 
Ausgangsdifferentialgleichung: 


dp < A 
at 2 Fl in—0 (Fr mit II<k=0).. . . (1.4). 


Man entwickle dann auch noch nach gan einer Stelle 9 = 9, und setze 9 — 9, = ydann erhält 
. man für (1.4) 


d? S o ee 9, 1 Sy e 
et Zu lh Hear ae erels beh): 


(Striche bedeuten hier und im Folgenden, daß die betreffende Größe nach y abzuleiten ist, eine 
darübergesetzte Null, daß sie an der Stelle 9 = 9, zu nehmen ist.) 

Es soll nun eine periodische Lösung der Kreisfrequenz zwei für die Differentialgleichung 
(1.5) konstruiert werden. Wir machen zunächst die Voraussetzung: 

a) Es sei in (1.5) Fo (g) nicht identisch Null. 

Zur Konstruktion der Lösung gehen wir mit dem Ansatz (1.3) in die Gl. (1.5)ein. In der 
Entwicklung (1.3) können nach (1.1) und (1.2), da rund « voneinander unabhängig sind, auch & 
und nals unabhängig voneinander angesehen werden. Man kann also verlangen, daß durch den 
Ansatz (1.3) unsere Differenzialgleichung identisch in & und n befriedigt wird. So erhält man 


dt 
durch Koeffizientenvergleich wegen were 41°C: 
k—1 k—2 
4% a, = Pan ı ; 10 R 
kı = ki, lij el Ri, let Bar - -» » . . (1.6) 
i=1 Av 


HI<i 
Rzı enthält nur ,;mit@ sk —3. 


!) Siegel, C. L.: „Über eine periodische Lösung im ebenen Dreikö blem.‘ 
an p g im ebenen Dreikörperproblem.‘“ Mathem. Nachr. 4 
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Im besonderen gilt für = 0 
— nıı 0=%-11 es +, — ME Da. 1-4 "45 14 +. . (1.7) 
und wenn man in der letzten Gl. (k + 1) anstatt % schreibt 
— Ei 4,0 = %,1 E = 0,51 3% Hi, a + Rerı,ı (1.8). 


Aus den beiden Gln. (1.6) und (1.8) lassen sich alle Koeffizienten a,;; rekursiv bestimmen: 
Glieder O-ter Ordnung treten nicht auf. Für k = 1 erhält man aus (1.6): 


122° 1 La) 
41 = qFPu 3 a, —_ et, 1 Fe LEE a ee (1.9) 
"und 
En = Es (99) RD EEE N en 9 (1.10). 


Die letztere Gleichung dient uns zur Bestimmung der ersten Näherung 9, für die Gleichgewichts- 
lage. Die Gleichung wird im allgemeinen nicht eindeutig lösbar sein, doch werden wir sehen, 


daß es für jede Wurzel von (1.10) unter der Bedingung Fio(p,) # 0 genau eine für hinreichend 


kleines & reguläre periodische Lösung (1.3) mit der erhaltenen Wurzel 9, als Anfangsglied gibt. 
Wir bestimmen die nächsten Koeffizienten: Für k—=2 (l # 0) erhält man aus (1.6) (beachte 
Id = 0 für 2 > k!) 


Ara, = N, +41 Fo+ 4, HH Fıt+ air fı,—- act > 5. (4,17). 
Aus (1.8) bestimme man a,, unter Beachtung von (1.9): 
o e) ons SCH [6 
— Foo = Fa+z Fur Fi, -ı +gFfu-1 Fu a (1.12) 


und 


Er digo — Fo+ Gy, Fi, 7 ao Fu nz 2a Ne nn. (ajo + 2 ayı 4.) | 


(1.13). 
+ 4yo (Fi: 41, Pr —1 ee] 
Schließlich schreiben wir noch hin 
41? ds, — B sl + IH Q2, Hr: a, Dt Ay; dl (1.14), 
u; 
> 2% Ay — Hut I 2% a4 IP @,-—ı+ nn Es ü3, = 
w et nn ae ll 


Er 2 Bei j-i' 41, ae > BIS i@o;'dı, | 


a] 
In dieser Weise werden aus (1.6) ke (1.8) in der Tat alle a;, rekursiv bestimmt, wobei stets 
(1.6) vor (1.8) anzuwenden ist. Offenbar ergibt sich so für jede Lösung 9, von (1.10) ein perio- 
discher Lösungszweig für die vorgelegte Differentialgleichung. Hängt die Differentialgleichung 
noch von einem 2. Parameter ß ab, so sind für «— 0 die Verzweigungsstellen durch 


Fio(90; B) = 0 
F (90, B) = 0 


gegeben. Setzt man hier weiter voraus Fo #0, so kann man auch an einer solchen Verzwei- 
gungsstelle die Frage nach der Existenz einer periodischen Lösung stellen. 

In diesem Falle liefert (1.12) eine zusätzliche Bedingung für die Koeffizienten 
nämlich 


) eo 0, ie Oo, - { 
M(g,) = Fa +7 Fa fı,+TF,m Fin Be RL ENCRO): 


Ist sie erfüllt, so erhält man weiter aus (1.13) unter Berücksichtigung von (1.9) und (1.11) nach 
kurzer Rechnung: 


1 or © | © o, 1 N 
ch hat ah + fu Fı-ı + HA) 


° 1 ) oO» je) Sur 1 oO Er 5 SV 72 
+Fy +7 Kaf, tan Fı)tz Fi -ı Faı +F,1'F%,-ı)+ Fi 


Z. angew. Math. Mech. 


108 Stellmacher, Über erzwungene nicht-lineare Schwingungen Bd. 34 Nr.3 März 1954 
en 2 a er nn 
ıl Or 

oder tm MHB=O 2.2.0: (1.17). 


Hierin ist M’= = die Ableitung der linken Seite von (1.16) nach @,. B enthält ebenfalls 
nur bekannte Größen. (1.17) ist unter der Bedingung Fj, 7 0 eine quadratische Gleichung zur 
Bestimmung von a. Wenn sie zwei reelle Wurzeln besitzt, so bestimmen diese beiden Wurzeln 
zwei periodische Lösungen, die von unserer Verzweigungsstelle ausgehen. Um die übrigen a; 
zu bestimmen, setzte man in (1.8), worin abermals k durch (k +1) zu ersetzen ist, aus (1.6) 
% 11,1 und @.+41,—ı ein. So erhält man 


la M’) = 2 it Dr (a1, aan.) (IB). 


(Ci, ®; sind bekannte Ausdrücke). Auf der rechten Seite stehen jetzt nur noch az; mit? 0 


und a;;, mit 7 <k. Ist a,, Wurzel von (1.17), so ist der Faktor von a0 in (1.18) bis auf den 


Faktor (Fjo)! die Diskriminante von (1.17). Ist diese von Null verschieden, so sind unsere 
Rekursionsformeln vollständig und damit alle Koeffizienten eindeutig festgelegt. Wir for- 
mulieren noch den Satz: 

Zu jeder Wurzel von (1.10) gibt es im Falle Fi, # 0 genau eine in der Umgebung von 
& —= O0 reguläre periodische Lösung von (1.4). Ist Fi, = 0 und (1.16) erfüllt, aber Fo 0, also 
4, eine doppelte Wurzel von (1.10), so gibt es genau zwei periodische Lösungen, die für « — 0 
zusammenfallen, wofern die Diskriminante der Quadratischen Gl. (1.17) negativ ist ?). Ist 
jedoch (1.16) an der Verzweigungsstelle nicht erfüllt, so gibt es dort keine in & reguläre peri- 
odische Lösung. 

Die Konvergenzbeweise holen wir im Abschnitt 2. nach. 


b) Fjo(g) ist identisch Null. 

Dieser Fall dürfte wohl das größere mechanische Interesse verdienen, wie wir später 
durch Beispiele belegen werden. Diesmal haben wir 9, zu berechnen aus (1.12), d. h. aus (1.16). 
Ähnlich wie oben im Falle der Verzweigung unter der Voraussetzung a) muß man, um die a;o 
zu berechnen, jeweils die a+ı,,1 aus (1.6) in (1.8), worin noch (% +1) für k zu schreiben ist, 
einsetzen. So erhält man zunächst entsprechend (1.17) 


; © ln Res o [e) j 
— te M(p)—= Fa +, KH harha) TS 


M’ ist die Ableitung der linken Seite von (1.16) nach 9,. (Dieser höchst plausible Faktor scheint 
aber nicht direkt ersichtlich, sondern ergibt sich nur durch unberirrtes Ausrechnen.) Allgemein 
erhält man entsprechend 


: 1 D, °, O, ) 
MM —z (a. Fi-4+ @,— Fı,4) +2 9%; 2, + 2a, ir F,5 01,7 | 
7 


+ Dan.) 

(D;0 enthalten wiederum nur Bekanntes). Damit sind die gewünschten Rekursionsformeln 

gewonnen. Die Eindeutigkeit der Entwirklung nach Vorgabe einer Wurzel von (1.16) ist 
gesichert, sofern der Faktor M’ von a;, in (1.20) von Null verschieden ist. 

Die Forderung, daß M’ 0 ist, ist also nichts anderes, als die Forderung, daß die Gl. 

(1.16) für 9 — 9, nach 9, auflösbar ist. Wie im Falle a) werden wir auch hier erwarten, daß 


120) 


man unter gewissen zusätzlichen Voraussetzungen an einer Verzweigungsstelle M'= 0 die 
Verzweigung der Lösung bis in die Verzweigungsstelle hinein verfolgen kann. Die zusätzliche 
Bedingung lautet diesmal nach Gl. (1.19) 


o ro o ) 6) 
No) =Fo+7 Fu Fa-ı4+ Fi Fr.) =U na rl 


während (1.20) mit M’—=0 und k—=2 nach Elimination aller ds, » nach einiger Rechnung 
ergibt: 


1 }; o Ir 220 ° o, o N 1 
5 % Me u fe +7 Fu Fe,—_ı+ Aa) +B, = zZ M"+ ao P+ B,=0 (1.23) 


wobei B, wiederum nur Bekanntes enthält. 


°) Ist auch #(9,) = 0 aber F’’(p,)£ 0, so erhält man entsprechend eine dreifache Verzweigungsstelle 
bei & = 0, nur falls außerdem (und außer (1.10)) noch M = 0 und B = 0 sowie M’ = 0 ist. Für den Fall MO 
siehe den nächsten Abschnitt b) 


ae re ra ee ee 
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Der Faktor von a}, erweist sich in der Tat als die zweite Ableitung der linken Seite von 
(1.16) nach g,. Danach erhält man analog (1.18) aus (1.20).(mit M'—=0Oundk>1) 


4,0° (Co M’’ + P) = Bekanntem. 


‚Wiederum sind im Falle des Nichtverschwindens des Faktors (a,, M’ + P), der sich wiederum 
im Wesentlichen als die Diskriminante der in a,, quadratischen Gl. (1.23) erweist, ae 
Rekursionsformeln vollständig. 


2. Konvergenzbeweis zu den vorangegangenen Entwieklungen 


Zunächst Fall a)mit fun #9. 
Der Konvergenzbeweis wird geführt im Anschluß an die in Anmerkung !) zitierte Sie - 
a Arbeit. Es ist zweckmäßig anstatt der Veränderlichen @ die neue Unbekannte 
= 9 — z 4: lie P— 9% — 9, einzuführen, sodaß also x mit Gliedern zweiter Ordnung 


ee Dureh Einsetzen in Gl. ’ 4 erhält man wieder eine Gleichung der Gestalt 


= tn) —0 


oder entwickelt: 


= = o 1 o, 
N 1 3 Flat ze Dfıkat-- Er (2.1). 


k=2 Kl 
Man setze dann 


Swen) r(2hıkut 2fai)=RiEn eu nn (22) 


Aus den Regularitätsverhältnissen der Funktion F und damit von f folgt die Majorantenab- 
schätzung 


a S 
R(®, |E|, B) *) M 3 D.= a @ ZN: 
BE (9 Zell) 3) 
wobei |&] = |n| = &, mit & bezeichnet ist, und ® für den Augenblick als unabhängige Variable 
behandelt ist. 

Man gewinnt aus den Rekursionsformeln (1.6) und (1.8) die Ungleichungen 


1 laul< lila EB DACH, ai NER 
il ES Ar Rat, 05 2 nen mo ee eh). 


Man multipliziere nun (2.4) mit a — L, und (2.5) mit |&go| = &, und summijere über alle k 
und l. So erhält man aus (2.4) und (2.3) 


&? 
101540 +M- + 


ron (A, passende Constante >0) . . . (2.6) 
und aus (2.5) und (2.3) 


b JJ Ye er Ber i ze Be 
fiel Bo 11 42 Pı+ ee, Bene; 2, le; = |mo||öro) 2 
Die letzte Beziehung wird weiter majorisiert durch formales Einsetzen von ®, aus (2.6) 

: 4 rd C+® C+(B +2) | 

AR U Eee EIER ji ME in ar), 
le (4: 2 T77(@ =) EN TE 


Addiert man nun (2.6) und mit passendem Faktor versehen (2.8), so erhält man nach Addition 
von & auf beiden Seiten und weiterer Abschätzung: 


2ER 
4(®8 +6)! Mi 0,(+ ©) EEE 
setzt man nun ® +& —=y so erhält man eine Beziehung 
aM SR DE REREURE 5 26 (2.10); 


Ofkı 


3 ° 
3) Hier ist = ha gesetzt usw. 


4) Das Majorantenzeichen $? bedeutet, daß die links stehende Potenzreihe in p kleinere Koeffizienten 
besitzt, als die rechtsstehende. M ist eine passende gewählte Konstante (hat nichts zu tun mit M in (1.16)). 
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Berechnet man jetzt aus der quadratischen Gleichung 4 Z(1—-,Z2)=M,;,( +2) 
die Reihe Z= u +£3...) (mit lauter positiven Gliedern!) 


so ist Z majorant zu y. Z konvergiert aber, solange die Diskriminante der quadratischen 
Gleichung Z2(M, +40,)—4Z+M,;t=0 D=(M, +4c,) ME —4<O ist, d.h. für 
4 


= M,(M;+4c,) 


Fallb)F=0. 
Wieder gehe man von der vereinfachten Form (2.1) aus. Diesmal ist es jedoch nötig, 
R etwas anders zu majorisieren: 


er BOE DRS 
R era M. En En RER TE TER 
Die Rekursionsformalen (1.6) und Eee vereinfachen sich nun zu 
AB ayı = af PER De ee (2.13), 
— M'’ a = Kto +— Sn Peek —ı Rrsı, ie 1%, a 4, —2) | 
(2.14) 
N 
sie N a nn 1, + Day, 4 | 
3,3 fi 


Hieraus gewinnt man die Ungleichungen 
4lay.| = Alla: +3 Bla-;t|Ril. . .... . - (2.15), 


I 1 7 er 
Map) S| Rural 4+G (Aal Br, | + fı, m RH.1)+ 23 0|%:|+23D |a-,.] (2.16), 
7 T 
wobei die 4A,, B;, C;, D; wohlbestimmte nicht-negative Konstante sind. Um nun eine Majo- 
rante für die Reihe ® = &, + ©, — 3 |azo| |&ro| + & a1] |&x | zu erhalten, multipliziere man 
1#0 


wieder (2.5) und (2.6) mit Mr — |£* und summiere über k und !. So erhält man aus (2.15) und 
(2.11) mit passendem A> 0 


48,1) A 2 EL UA ER 
NACPLEDHL I I oT 


Um entsprechend mit (2.16) verfahren zu können, beachte man, daß in (2.16) stets % > 1 bleibt, 
danach benutze man Majorisierungen 


Silkaltr ge: Zilk; Zlhlliuitt 


2.1. 


k=2 er k= ze er 
Damit erhält man schließlich aus (2. “ unter Benutzung von (2.11) 
Me 2® +:08 N 
M'4: D,1} - 1—4£ (4: 12 | M, an nn Crema er (2.18). 
Sr neue Majorante hierfür erhält man unter Benutzung von (2.17), das rechts formal in (2. 18) 
für ©, eingesetzt werde: 
M, 5 sd +L® 
M':4® 2 Aldc+ or M, AAN 
Taerar 4m. ++ ta li aldr 20) 


addiert man hierzu nach Division durch M’ (2.17) unter Berücksichtigung von (2.11), so erhält 
man 


1 
104 Me Ä ana. 
en: NT 
er : 1 
B+o+@, (@+ ro o) eh, 
17 i MM; S 
—(®+I)% ı 


1- (+2 2). 
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Addiert man auf beiden Seiten 4 £? und setzt man 


1 
so ergibt sich 
Sry 
a END Ann rule rel 


Diese Abschätzung ist wieder von der Gestalt (2.10) wie dort schließt man auf die Konvergenz 


für &d< —— ———— , 
M, (M;+4 0,) 

Erineinfaches-Beispiel 

Um zu erkennen, daß in einfachen Fällen bei derartigen Störungsentwicklungen der Kon- 
vergenzradius durchaus nicht besonders klein sein muß, was für numerische Berechnungen ja 
sehr wichtig ist, schätzen wir in einem besonders einfachen Falle erzwungener Schwingungen 
den Konvergenzradius explicite ab. Zugleich erkennen wir dabei, wie geeignet die Siegel- 
sche Abschätzungstechnik auch für derartige Abschätzungen ist. Sei vorgelegt die Differen- 
tialgleichung 


d’p 2 2 
18 +n(9+bp) +ecoswot=0 
oder mit den Bezeichnungen 


2 
n C & A 
0 —4; dt. — &; Bene ne 


oat=T —,=& —, = 
E) w? 3 2 


£ - d? @ 
Eu et! ak? Fa + &10 (9 +5?) hy Stra). 
Der Ansatz (1.3) liefert hier 
7 
4,5; Er 3,.2=0; med wel m—0... 


und allgemein für k>1 
R 1 = %ı,ı Farb-(ar 2,11 + &-2,142) + 
—.4,—0 = b = (H—ı-#— ee) a 
‚ 
Geht man nun hier wie im allgemeinen Falle des Konvergenzbeweises für den Fall a) vor, so 
erhält man entsprechend den Abschätzungen (2.6) und (2.7) 
Bl ac+t(d, +58) 
8, b@: 
8 +Dd=Blac+Bc+@b(CcH1), 
wie bei (2.10) genügt eine Majorante Z von ® der Gleichung 
Ze alöE -ZE + ZUIbIE 1, 0 See AN 
und konvergiert solange die Diskriminante dieser in Z quadratischen Gleichung <0 ist, d.h. 
für 
BE RE naar Henn (BD) 
Sei nun & = x die absolut kleinste Wurzel von (B), nachdem darin statt < das Gleichheits- 
zeichen gesetzt ist. Dann ist die Majorante Z und damit unsere Reihenentwicklung konvergent 
für&<|®]. Man errechnet, 2 |a b| = gesetzt, 


el 1 »\l1+42)) 


und sieht daran, daß für kleines A = 2|« b| unsere Entwicklung sogar bis dicht an die Resonanz- 
stelle &—= 1 konvergiert. Allerdings nimmt der Konvergenzradius bzw. die Schranke für den 
Konvergenzradius, die wir somit gefunden haben, mit wachsendem A monoton bis auf Null ab. 


n? 
Andererseits weiß man, daß für negatives b oberhalb der Resonanzstelle (a: rhurg — Pr — ı) 


eine Verzweigung unserer periodischen Lösung existiert, die mit @b— 0 sich der Resonanz- 
stelle & —= 1 nähert. Daher kann unsere Abschätzung des Konvergenzradius in diesem Falle als 
befriedigend bezeichnet werden. 
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3. Stabilität der gewonnenen Lösung 


Um die Stabilität der gewonnenen periodischen Lösung p(z) unserer Bewegungsgleichung 
zu beurteilen, fügen wir, grob gesprochen, eine kleine Störung zur Lösung hinzu und prüfen, ob 
diese Störung klein bleibt oder nicht. Zu dem Zwecke variieren wir die Lösung o(r) d.h. wir 
denken sie eingebettet in eine einparametrige Schar von (im allgemeinen nicht periodischen) 
Lösungen B(t, «, e) von (1.5), die als regulär analytisch in a, vund & vorausgesetzt werden. Es 
sei lim B(r, &, e) —=g(r). Durch Differenzieren von (1.5) in das die Lösungsschar B(t, a, &) ein- 

e>0 


gesetzt zu denken ist, nach dem Scharparameter & erhält man, da nur p von & abhängen soll, 


oB BR 
(vonn = u(t, &) gesetzt wird) 


d2 { oF 0, Ir 1 zv 
Zr +Ae)a=0; =, = D Furtat+ 2 Frl ytzy Pu; | (3.1). 
ee ED RE ET h(t) = I hrı Sn gesetzt. | 


(3.1) ist die Jaco bische Derivierte der Differentialgleichung (1.5). 

Wie üblich, nennen wir die in den vorigen Abschnitten gewonnene periodische Lösung 
von (1.5) instabil, falls es mindestens eine Lösung der zugehörten Derivierten (3.1) gibt, die, 
absolut genommen, für hinreichend großes 7 beliebig große Werte annimmt. Umgekehrt nennen 
wir die periodische Lösung stabil, wenn alle Lösungen der dazugehörigen Derivierten (3.1) die 
Eigenschaft besitzen, daß der absolute Betrag der Lösung u des Cauch yschen Anfangswert- 

du 
(= 
gegebene Schranke ö herabgedrückt werden kann, wenn nur c} + c$ hinreichend klein gewählt 
wird. Um das Kriterium anzuwenden, sucht man Lösungen von der Gestalt 


N == eiut.ny s Me > DI ee at re a ee er (3.2), 


problems von (3.1) für die Anfangswerte u(0) = c,, — (, für alle Zeiten unter jede vor- 


wobei wiederum v periodisch von einer ganzzahligen Periode sein soll. Die Durchführung dieser 
Aufgabe ist wohlbekannt und führt zur Bestimmung von u auf die bekannte Hillsche unend- 
liche Determinante. Ist u reell, so herrscht Stabilität im obigen Sinne, ist u rein imaginär, so 
herrscht Instabilität (R(r) in (3.1) sei reell!). 

Wir wollen hier einen von der Hillschen Methode abweichenden Weg einschlagen, da 
wir ja die Lösung p und damit auch die Funktion Ah(r) aus (3.1) als Potenzreihe in &, n gegeben 
haben. Will man nach der Hillschen Methode die Glieder niederster Ordnung einigermaßen 
genau bestimmen, so hat man gleichwohl in der Hillschen Determinante schon recht viele 
Glieder zunehmen. Auch ist es für uns nur folgerichtig, weiter an der Potenzreihenentwicklung 
festzuhalten. Das wird uns die Berechnung in der Tat vereinfachen. Wir vermeiden also den 
Begriff der unendlichen Determinante und stellen ähnlich wie im Vorangegangenen einfache 
und für praktische Berechnungen brauchbare Rekursionsformeln auf. Auch der Konvergenz- 
beweis ergibt sich zwanglos. 

Geht man mit dem Ansatz (3.2) in die Differentialgleichung (3.1) ein, so erhält man die 
neue Differentialgleichung 

d?v 


.. dv ; R ; 
eian 2i HZ +Fvchle) =) U Rn ee 
Hierin soll nun 
I PYT7 ak = 2 ur Ge ee (3.4a) 


so bestimmt werden, daß die Differentialgleichung (3.3) eine mit der Kreisfrequenz zwei peri- 
odische Lösung besitzt. Wegen der Linearität und Homogenität von (3.3) ist vin jedem Fall 
nur bis auf einen zeitunabhängigen Faktor bestimmt, der aber natürlich noch von « abhängen 
darf. Wir wollen annehmen, was sich bald bestätigen wird, daß der zeitunabhängige Summand 
Uno &;0 der Entwicklung von v von Null verschieden ist, so daß wir ihn zu Eins normieren 


können. Demgemäß machen wir den Ansatz 


v=1l+ 2 ud Chin) vn ei, 


ie ee 
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Wir setzen voraus, daß Ag = Fü = 0 und ho = Fi = (ist 5). Das bedeutet, daß wir uns nur 
auf sogenannte ganzzahlige Resonanz beschränken und Resonanzfragen höherer Ordnung außer 
Acht lassen, wie es für den hier interessierenden Bereich von Schwingungen hoher Erregerfre- 
quenz nur natürlich ist. Geht man mit dem Ansatz (3.4a) und (3.4b) in (3.3) ein, und verlangt 
wieder identisches Verschwinden in & und n so folgt zunächst bei Betrachtung der Glieder 
Null-ter Ordnung u, —= 0. Weiter durch Annullieren des Faktors von &,; für 1 #0 


4, = —Al 93 Yo t ut ur er — Svuimzıky 0: - - (3.5). 
i 7 ’ ij 
Entsprechend erhält man für } = 0 die Bestimmungsgleichungen für die u; 
k—2 


2 ı m = 2 iii + rot I ur er. 


Hierin setze man noch &k +1 für k 


k—1 
2 Kı Ur = - Wi; Uk-+1—i + hr-ı, ) + DE UR-1—k', Ir Apr ee ee (3.6). 
ıi= 
Im einzelnen erhält man dann für die Glieder erster Ordnung aus (3.5) unmittelbar 
1 De cr 
AN zil — hrs = Fi 2; Ve )) De ie ee (3.7) 
ferner für k==2, 1=0 

j u“ ee ho + Aı v1, —1 -- hi, | ’ı1 BE a har a SE a ee ae (3.8) 


d.h. bei Beachtung von (3.7), (3.1) und (1.9): 
2, 1 En - en > 1 Sn Dir ! 
= Fn+z Fiıfı,-t F1.-1-Fu) + Fu Fi, = Mg) ee, (29). 


Hier tritt wieder die Empfindlichkeitsgröße aus (1.19) auf. Das steht mjt unserer Erkenntnis 
im Einklang, daß die Bedingung E = 0 für «— 0 Verzweigungsstelle für die periodischen Lö- 
sungen von (1.5)ist. Denn wenn eine periodische Lösung von (1.4) eine Verzweigungsstelle be- 
sitzt, so hat bekanntlich die dazugehörige Jacobische Derivierte (3.1) immer eine peri- 
odische Lösung. Doch gilt in unserem speziellen Falle darüber hinaus: Wenn für «— 0 eine 
periodische Lösung der Jaco bischen Derivierten existiert, so besitzt die Ausgangsdifferen- 
zialgleichung dort immer eine Verzweigungsstelle, wie aus (3.9) unmittelbar abzulesen ist. 


Wir können für jede periodische Lösung von (1.4) feststellen, ob sie stabil ist oder nicht, 
indem wir in der angegebenen Weise u berechnen und prüfen, ob diese Größe positiv oder nega- 
tiv ist. Das Verfahren versagt natürlich an der Verzweigungsstelle selbst, da dort die peri- 
odische Lösung nicht mehr konvergiert. 


Nachdem u, durch (3.9) festgelegt ist, kann man rekursiv alle Koeffizienten u;, yı ein- 
deutig festlegen aus (3.5) und (3.6). Die Glieder niederster Ordnung seien noch angegeben 


4m = — Alta + har + vu Ayımı + dı,-1 hy,ır | 
2 1 a = Rat Ya hr, 4 9%, hut thai tt,ihı | ee 3i10) 


j; 
—= he 5 Aıkı, + 6% (Ag, —ı Aıı + Aaı hı, -ı) 


o 

5) Im Falle a) des Abschnitt 1 Fl, = hy 0 ; hoo = 0 führt unser Ansatz (3.4b) nicht ohne weiteres’ zum 
Ziel. In diesem Falle setze man &* = &; n7*? = n und entwickle nach &*, 7* dann ist in unmittelbar verständlicher 
Bezeichnungsweise wieder hf, = 0; h%, 0. Danach verfahre man wie oben. Auf diese Weise erhält man die so- 
genannten halbzahligen Lösungen (Kreisfrequenz 1). Ist auch Ayo 0 (ein Fall, der aus dem Rahmen der vor- 
liegenden Arbeit herausfällt, aber mühelos miterledigt werden kann,) so entwickle man zweckmäßig wieder nach 
&*, n*. Dann hat man zunächst durch Koeffizientenvergleich von &#: u3 = hd = Iive. Danach geht alles wie oben. 
Die Verzweigungsstellen bzw. Stabilitätsgrenzen im Grenzfall «— 0 erhält man diemal in der Form hf, = u = 1? 
d.h. YA, muß eine ganze Zahl sein. Dann ist 4 nach (3.2) periodisch. Anstatt (3.5) erhält man dann 4 (1 + wo) Urı 
= — 412 m-iım.... 

6) Ganz entspr. wie in diesem Abschnitt kann man vorgehen, um den Satz zu beweisen: „An jeder 
Stelle der Stabilitätsgrenze liegt (falls dort die vorgelegte Bewegungsgl. regulär ist) eine Verzweigungs- 
stelle. Liegt eine Verzweigung in eine gerade Anzahl von Ästen vor, so ist umgekehrt auch die Verzwei- 
gungsstelle Stabilitätsgrenze“. — Zum Beweis hat man nur statt für a = (0 nahe an der Verzweigungsstelle 
zu entwickeln. 


8 
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4. Konvergenzbeweis zur vorangegangenen Entwicklung 
Um den Konvergenzbeweis zu führen, ist es zweckmäßig, die Gl. (3.3) etwas umzuformen. 


Zunächst setze man 
seltltnluotn astra . (4.1), 
mit imei-n tr m re We 
Damit erhält man aus (3.3) 
> 5 ed [Pr m I ed +(1+v)2» p+P)+w(hr) - 7 — Pp (2, +P)) (4.3), 
+htr)—=0 


wo nun w und p als Unbekannte », und v, als bekannt angesehen werden können. Dabei ist 
gesetzt 
dv, 


h(r) = AS, en +(d + v1) (Ar) — 97) IE RR EINEN (4.3a). 


Wieder suchen wir eine Majorante für die w, p entsprechenden Reihen 


Wie oben ist &,, = |&*t! Be = o% = L, zu setzen. 
Da h(r) regulär analytisch in £, n vorausgesetzt wird, so gilt dasselbe auch für ho) daher 
BD u, ee ee 
1 ch 
Läßt man hier links die Glieder erster und zweiter Ordnung fort, so erhält man für die übrig- 
bleibende Funktion 
= (6 —: 
H, PA a an ee 


Aus den ln. (3.5) und (3.6) gewinnt man ha folgenden Ungleichungen. 


4 \vg:] 245 z en ii + |Arıl + &. 1 | Mar) — 


> mir ‚ıllasuel;‘ & 22 Tess 
2 al: < a |; ka 


Multipliziert man (4.7) mit |&.,| = &, und (4.8) mit &;, = &, und summiert über %, 1, so erhält 
man unter Berücksichtigung von (4.3), (4.4), (4.5) und (4.6) 


Aw KmlE+P)W+AP Kol + )IE + +Llaal+ le, all) 22 | P+ P2) 
+W(H + nP + Pf2w|E+P)) +H 


M 1 ö 
N u TABU TP)+P@e + PO)+ Pe Fa" + u: 


Wu er pre (PM BR 
ERST To; Sue PT ZEN; 
Die Richtigkeit der letzten Majorisierung erkennt man, wenn man rechts den Zähler entwickelt, 


und vom Nenner zunächst nur 1— c,& berücksichtigt. Danach ist nur noch festzustellen, daß 
auch ein Glied von der Form WP?rechts enthalten ist. 


Entsprechend erhält man aus (4.8) 


(4.9). 


M ex c3 ME 

I) P2 3 ; 

2 ml PP) PH 1 Zt Fe 
P a WM Men 
N Ä &0 
ale Iren 


Setzt man hier für W/£ die letzte in 


| ier fü (4.9) gewonnene Majorisierung ein, so sieht man die Richtig- 
keit der Majorisierung 


By, mE HM + BR +18 
ST GTZ, (WeErBE 
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Schließlich durch Addition von (4.9) hierzu 


ie: Se er Hear e): 
A en ir Lt u Re 
er 


Setzt man nun = + * -+-£so hat man 
| ty? 
4 M a ee NE ROHR 
vu Mk 2 (4.11). 
Also wieder die Form (2.12) und schließt wie dort auf Konvergenz, wenn nur 


£ 4 
= M, (M&+ 6) i 


5. Beispiel: Pendel mit oszillierendem Aufhängepunkt, lineare harmonische Erregung 


Die gewonnenen Ergebnisse sollen zunächst angewandt werden auf die Bewegung eines 
schweren Pendels, dessen Aufhängepunkt geradlinig und periodisch beschleunigt wird. 


2 
98 + mgssinp+ ms Bl)sin @—)=0. ET ME EL OLE: 


Es bedeutet (s. Bild1) 9 den Winkel zwischen 

Pendel und Schwererichtung g(lJ|=g), ö dnw 
spitzen Winkel zwischen Pendel und Beschleuni- a, i 
gungsrichtung B(|®] = B(t)); ©, m, s, Isind Träg- 
heitsmoment, Masse, Schwerpunktsabstand und 

reduzierte Pendellänge des Pendels. Die Erschüt- 
terungsfunktion B(t) sei zunächst harmonisch: 


Di) 0.4.6005 onen. 293 102), 


Um die Gl. (5.1) dimensionslos zu machen, setze 


man 
| ng g 
a a 


und erhält 


2 
= +o®sinp +4%cos 2 7-sin (pP — ö) gi, 
dr (5.2b). 


Bild 1. Pendel mit geradlinig beschleunigtem 
Es kommt nun darauf an, auch den Koeffizien- air 
ten 4 A/l des Beschleunigungsgliedes von a ab- 


hängig zu machen. Um den Anschluß an die vorausgegangenen Entwicklungen zu gewinnen, 


setzen wir 

Pe —2ob; = | 

| Yig 
Danach wird b als von & unabhängig angesehen. Eine solche Zuordnung ist erlaubt, weil wir 
die Gültigkeit der vorstehenden Reihenentwicklungen für ein festes Intervall |x| < a, nachge- 
wiesen haben. (5.3) ist der einfachste Ansatz, der die Möglichkeit gibt, für 9, allgemeine Werte 
zu erhalten. Der Ansatz ‚‚A/l von a unabhängig‘ liefert nämlich sofort g, = ö (Pendel in der 
Erschütterungsrichtung) und der Ansatz A/l = a?-b* führt auf 9, —0. Wir wollen ja aber 
gerade solche Bewegungsvorgänge beschreiben, bei denen die mittlere Auslenkung aus ‚der 
Gleichgewichtslage möglichst stark ist. Da alle unsere Entwicklungen vorzugsweise für kleines 
«& gültig sind, besagt (5.3) physikalisch, daß wir auch das Verhältnis A/lals einigermaßen klein 
ansehen. 
In der Bezeichnungsweise der Gl. (1.4) haben wir jetzt anstatt (5.2b) 


. (5.3). 


2 . 
U REERR? ee ch D.2E), 


F„=sng, Fu=Ff,-=bsin(,— 6), übrige Fu =0. 
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Als Mittellagengleichung des Pendels ergibt sich nun aus (1.16) in „‚Null-ter Näherung“ 
2 


My) =sin + in2(—d)=0. u he re ae 


Das ist dieselbe Bestimmungsgleichung wie sie für unser Pendel durch. andere Näherungs- 
rechnungen wohlbekannt ist, und von verschiedenen Seiten r) überein ine bestätigt 


4A? w? 

wurde. Aus (5.5) liest man ab: Für sehr kleines b? (d.h. nach (5.3) für kleines G 
zwei Gleichgewichtslagen nahe bei den Gleichgewichtslagen des ruhenden Pendels. Bei wach- 
sender Erregerfrequenz kommt noch ein weiteres Paar von Wurzeln der Gl. (5.5) hinzu, so daß 
für großes b? vier Gleichgewichtslagen vorhanden sind, nahe der positiven und negativen Er- 
schütterungsrichtung und der beiden dazu senkrechten Richtungen. 

Auch die sich hier ergebende erste Näherung für den Bewegungszustand des Pendels ist 
wohlbekannt und wird hier nur bestätigt. Man erhält nach (1.9) und (5.4) 


hat man 


P= 9 +F sin (9) cos 2r a 


Eine Gleichung, die übrigens besagt, daß das Pendel für kleines A/l jeweils infinitesimale 
Drehungen um den ‚‚Stoßpunkt“ des Pendels ausführt, (Der Stoßpunkt besitzt bekanntlich 
den Abstand } vom Aufhängepunkt des Pendels). Man beweist dies leicht, indem man zeigt, 
daß für kleines A/l die Kurve, die der Stoßpunkt in einem Inertialsystem beschreibt, mit ihrer 
Tangente im Stoßpunkt auf den Aufhängepunkt des Pendels zeigt. Die Bewegungskurve des 
Stoßpunktes im Inertialsystem wird während einer Schwingung von der Mittellinie des Pendels 
eingehüllt (wofern die Bewegung des Pendels durch die Näherung (5.6) beschrieben wird). 

Zur Beurteilung der Stabilität haben wir nach (3.9) die Empfindlichkeitsgröße u oder 
näherungsweise nach (3.9) für kleines 1a =& YM’ auszurechnen. Man erhält etwa durch 
Differenzieren von (5.5) oder durch direktes Ausrechnen von (3.9) 

2 


M = cos +2 00820). Er N 


Wenn diese Größe positiv ist, so ist die Bewegung für hinreichend kleines «a stabil, für M’<0 
labil. Betrachtet man in (5.5) b?als durch diese Gleichung gegebene Funktion b? — f(9,, 6) und 
zeichnet man die Höhenlinien dieser Funktion in der b, ö-Ebene, so erhält man eine Schar von 
Sinuslinien, deren Enveloppe nach (5.7) die Stabilitätsgrenze für &— 0 ist, weil (5.7) die Ab- 
leitung von (5.5) nach 9, darstellt. 

Durch (5.7) ist nun in der Tat die Stabilität in der Weise festgelegt, wie es einige Autoren 
auf Grund direkter mechanischer Überlegungen gefunden haben ®). Jedoch erhielten später 
Klotter u. Kotowski (l.c.°) ein abweichendes Ergebnis auf Grund einer Inte- 
gration der Bewegungsgleichung, die dabei jedoch in unzulässiger Weise linearisiert wurde, so 
daß das von Klotteru. Kotowski (l.c. S. 296) in Gl. (31) besonders erwähnte Glied 
verloren ging. — Leider sind auch die Stabilitätsbetrachtungen von Haacke?°), die wesent- 
lichan Klotteru.Kotowskianknüpfen, aus dem nämlichen Grunde einer nicht konse- 
quenten Linearisierung der Bewegungsgleichung nicht stichhaltig. Anstatt der Matthieu- 
schen Differentialgleichung (bei Klotter u. Kotowski Gl. (4b)) ist eben die unserer 
Gl. (3.1) entsprechende Hillsche Differentialgleichung zu diskutieren 1), 

Grundsätzlich scheinen die Stabilitätsverhältnisse verfälscht zu werden, wenn man den 
Schwinger von vornherein als rheolinearen Schwinger behandelt, und zwar um so mehr, je 
stärkere Veränderungen der Gleichgewichtslage des Schwingers durch die Erschütterungen in 
Rede stehen. 

Genau so einfach läßt sich der Fall erledigen, in dem der Aufhängepunkt des Pendels eine 
elliptisch-polarisierte harmonische Bewegungausführt. Zerlegt man 
die auf das Pendel wirkenden Beschleunigungskräfte, nach Komponenten in Richtung der 
beiden Hauptachsen derjenigen Ellipse, auf der der Aufhängepunkt des Pendels geführt werde, 


?) Eine kritische Sichtung älterer Ergebnisse findet man bei Klotter und Kotowski A 
Math. Mech. 19 (1939) 8. 289—296. Dort findet man auch weitere Literaturangaben. 
®)G. Greenhill: Rep. and Mem..Nr. 238, 1915; P. Hirsch: Z, angew. Math. Mech. Bd. 10 
(1930), 8.41—52; K.Klotter: Jahrbuch d. Deutsch. Luftfahrtforschg. 1933, S. III, 3/11. 
9) W. Haa cke- Z.angew. Math. Mech. 31 (1951), 8. 161—169 u. 8. 333 —338. 
10) Die von Klotter u. Kotowski angegebene Matthieusche Differentialgl. ist mit unserer 
(3.1) beim vorliegenden Beispiel nur im Falle ö = 0 identisch, weil in diesem Falle %=0 und, =n exakte 


Lösungen von (5.2) sind Aber natürlich ist die Lage dadurch festgelegt. Einen indifferenten Bereich wi i 
Haacke, I. Mitteilung gibt es nicht, N re 


Z. angew. 


ee  Stellmacher, Über erzwungene nicht-lineare Schwingungen 117 


so gewinnt man als Bewegungsgleichung des Pendels 
2 
= +mgssing+ms A, w!cos wtsin(p — d)+ms A,w®sinwtcos (—d)—=0 (5.8). 


Dimensionslos geschrieben mit den Bezeichnungen und Festsetzungen (5.2a) und außerdem 
A 


4 —2bo; u N (24-47,7 zen ee 3 . 13:9) 
1 l YI:g VI-g 
gesetzt, erhält man die für unsere Entwicklungen passende Form 
de 
zu + Fakt Farlu+ Pass Se Ver DEE RERIOFLDN 


Diesmal berechnen sich die Koeffizienten F},, wenn man noch setzt 9, = b,sin (9 — Ö); 
9 — b, cos (P — 8), i = Y—1, 


1 1 
F,=sin 9; RE H,-a4=n-9 un: 
Jetzt erhält man zur Berechnung der Gleichgewichtslage wieder nach (1.16) 
1 | 
M(9,) = sin 9 -+ 5 (& — ba)sın 20 0) = 0 ea 


(durch P. Hirsch zuerst angegeben und diskutiert). Der Stabilitätskoeffizient berechnet 
sich wiederum einfach durch Ableitung nach 9, 


‚ 1 E 
Mg) = cos +7 (di — b2) cos 2 (m — 8) - ma --(3-.13). 


Im Falle b, = 0 erhält man natürlich wieder die entsprechenden Formeln für das Pendel, dessen 
Aufhängepunkt linear harmonisch bewegt wird. — Für kleines « bedeutet wiederum M’>0 
Stabilität und M’< 0 Labilität. (Dieser einfache Zusammenhang scheint noch nicht bemerkt 
worden zu sein.) 

Um Korrekturen für die Gleichgewichtslagen zu bestimmen, hat man nach (1.20) 
und (1.21) die a;, zu berechnen. Nach Gl. (1.20) und (5.11) ergibt sich a, —= 0. Für die in « 
quadratische Korrektur der Gleichgewichtslage «? a,, erhält man im Falle der Gleichungsform 
(5.10) nach einiger Rechnung 


1 IR 
iM | (Fi: ‚a ie Ab ı) (Bat Mırkı —ıgFu‘ hm) 


+ Pu: PR, | 


Win 2m tim) N... 010 


— 


15 
m 
[ep) 
S 
—— 
I 


+ b3 sin? (9, — Ö)) — 


28 


16 (bi sin? (9, — 6) + b} cos? (9, — 2 


— b? sin? (9, — 6) + b2 cos? (9, — 6) sin 
Im Falle b, — 0 gilt die letzte Gleichung wieder für lineare harmonische Erregung. 

Interessanter ist jedoch die Korrektur für die Stabilitätsgrenze u, die von erster Ordnung 
ist. Sie zu berechnen benutzen wir (3.10), (3.1) und (5.11) und erhalten zunächst 

a, m hzı =. ha, 1 = 0 
dann folgt aus (3.10) und (3.1) sowie (5.11) 
Det er | e 
U = — Fu Fi, Ser ra cos? (9 — 6) +b3sin? (my —6)) - - - (5.15). 

Wie man sieht, liefert diese Korrektur stets eine Stabilitätsverschlechterung. 


Nach (3.2) kann man u = aw-+0?w-+ '' deuten als die reduzierte Kreisfrequenz der- 
jenigen Schwingung, die nach einer kleinen Störung der rein periodischen Bewegung des Pen- 
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_ dels überlagert wird. Nach (5.2a), (5.7) und (5.15) errechnet man als Frequenz dieser langsamen 
Schwingung im Falle b,—=0 


wo _® RN ER rn 
es = 4„(emt+% k)=5, rm (9 — 9) SEE Po RR 


TE gT 


so erhält man für nicht zu kleines A/leine Korrektur, die man an einer passenden Apparatur 
nachprüfen kann. “ 


Als weitere Möglichkeit für die Anwendung unserer mathematischen Entwicklungen 
bietet sich der Fall dar, daß ein schweres Pendel in allgemeinerer Weise periodisch erregt wird, 
z. B. derart, daß der Aufhängepunkt des Pendels durch ein Getriebe geführt wird, das seiner- 
seits durch einen Antriebsmotor bewegt wird. Das sei hier noch kurz für den Fall erläutert, daß 
der Aufhängepunkt des Pendels durch ein Schubkurbelgetriebe geradlinig periodisch bewegt 
wird. Mit den Bezeichnungen wie sie aus Bild 2ersichtlich sind, erhält man aus L= 22 +r1r — 


— 2xrcoswt ” 


TE TEN. 
= r (cos ot) +YJR — r?sin®wt=rcoswt+ »yı-(2) sinio tra 
und durch Entwicklung der Wurzel 
18% Lie‘? ar 
A Be ie 6 BURN, 5.18). 
s=ltr|onos Du A ot (2) sintwt 1elz) sine ot ( ) 


Hieraus endlich durch Fourier entwicklung und Differenzieren: 


2 4 6 
za relz)osort2[z) o2ott zz) [eos 2 ot cos don (7) ig! 
519 


Setzt man dies für B(t)in (5.1) ein, so ist es naheliegend, in Analogie zu (5.3) diesmal r/L und »/l 
als kleine Größen anzusehen und demgemäß als Ergänzung zu (5.2a) den Ansatz 


zu machen, was im allgemeinen auch physikalisch 
sinnvoll sein dürfte. Damit haben wir uns wieder 
unserem allgemeinen Ansatz (1.4) angepaßt. 

Die Korrekturen, wie sie für die Abweichung 
einer solchen Bewegung des Aufhängepunktes gegen- 
über einer harmonischen Bewegung des Aufhänge- 
punktes anzubringen wären, wirken sich lediglich für 
die Korrektur (5.14) der Gleichgewichtslage aus. Die 
Korrektur a,, bleibt Null und auch die Korrektur u, 
der Stabilitätsgröße bleibt unverändert. 

Alle diese Ergebnisse über das erschütterte Pen- 
del kann man naturgemäß gewinnen, ausgehend nur 
von einem „Duffingschen Ansatz“ der Gestalt 
(3.6) für die periodische Lösung, deren Stabilität dann 
weiter untersucht wird. Es dürfte jedoch rechnerische 
Schwierigkeiten bereiten, etwa durch Hinzunahme 
mehrerer Cosinusglieder höherer Ordnung diesen An- 
satz zu verfeinern, um Korrekturen entsprechend 


Bild 2. Schwerependel, dessen Aufhängepunkt durch ’ 1 
oin Botmunlaurbeiseiriohessrane "en den obigen zu gewinnen. Nahezu ausgeschlossen darf. 


“es erscheinen, gleichzeitig ein rechnerisch brauch- 


bares Verfahren auf einen solchen Duffingschen Ansatz zu gründen und dessen Kon- 
vergenz nachzuweisen. 


Die in der vorliegenden Note gegebenen Konvergenzabschätzungen dürften sich zum 
Teil dadurch umgehen lassen, daß man durch passende Umformungen den Anschluß an 
die bekannten Entwicklungen von Poincare herstellt, und sich dann auf dessen Abschätzun- 


Y 3 ET, j 


. Se hwir [Q 


. ich a a Eich dieser Fan ist inw unserem 2 
it enthalten, so daß auch dafür hier alle Formeln bereitgestellt sind. 


/ Es dürfte sich lohnen und keine Schwierigkeiten bereiten, die vorliegende Ontersuieh 3 
auf mehrere Freiheitsgrade auszudehnen, um somit etwa den Fehler eines Kreiselkompasses 
auf schlingerndem Schiff und ähnliche Fragestellungen durch direktes Zurückgreifen auf die 
en Bewegungsgleichungen zu behandeln. 


er Klotter, Forschg. a. d. Geb. d. Ing.-Wes. Bd. 12 (A) 1941, S. 214. 


Eingegangen am 16. Februar 1953. 


KLEINE MITTEILUNGEN 


en sukunen und Bernoullische Zahlen 
Die moderne systematische Darstellung !) der Potenz- 


summation erklärt zuerst die Bernoullischen 


Zahlen B, durch B,=1 und die symbolische „. Re- 
kursionsformel 


Baer bei", .W, 
dann die Bernoullischen Polynome B,(x) durch 


Sichtlich sind links n und & vertauschbar, also gilt 
Sr (n) + (S(a) + n)E = (Sr(@) + (Sm) +a)k (6) 


Dies fassen wir als eine Identität von Polynomen in « 
auf und vergleichen die Koeffizienten der in x linearen 
Glieder. Bezeichnen wir dafür das lineare Glied von 
S;(x) mit B;- x, so ergibt sich, wenn wir den Koefli- 
zienten des linearen Gliedes links sofort wieder sym- 


4 die symbolische ?) Definition bolisch schreiben: 
= Be) kBäs seen er. Ban Be + (1) Se, 
4 Schließlich gewinnt man für die Summe womit zunächst (4) (mit k — 1 statt k) erreicht ist. Es 
Fr Syn) =0k +1 +--:+(m—1)k (k>0) . (3) macht aber keine Mühe, durch » = 1 auch (1) daraus 
zu entnehmen. Nach (3) ist nämlich 

der k-ten Potenzen der ersten n ganzen Zahlen die ale . 
Formel Su VEaT arbeite 2 re) 
Fit S Bra WM) By 4 Die Definition (2) der Bernoullischen Polynome 
»(n) = De Te, - (4). wird durch (7) nahe gelegt. Dieser Weg zur Potenz- 
= summation hat den großen Vorzug, zwangsläufig aus 


Diese Darstellung ist vom heuristischen Standpunkt 
unbefriedigend, weil man darin kaum noch erkennt, 
wie das Problem der Potenzsummation zur Entdeckung 
der Bernoullischen Zahlen zwingt. 

Aus diesem Grund ist vielleicht folgender Weg zur 
Auffindung der Potenzsummenformeln von Interesse. 
Es sei schon bekannt, z. B. aus der Theorie der arith- 
metischen Reihen, daß S;(n) ein Polynom in » ist. 
Dann bilde man mit ganzem positivem n und x: 


dem binomischen Lehrsatz und der Definition der 
Potenzsummen zu den symbolischen Formeln zu führen. 

Auch der folgende zweite Weg scheint mir der Be- 
achtung wert. Setzen wir S;(n) wieder als Polynom 
voraus, so können wir es nach n differenzieren und er- 
halten für die Differenz der Ableitung: 


AS; (n) = Sm +1) — 80) = Z.{Sxn +1) - Sul) 


d 
Sn +)=Rk+l +-.-+(m— 1% a rinnen (kl)e 
n Glieder (5). Dieselbe Differenz hat aber auch das k-fache der (k —1)- 
Ink +tH(n +1) +. +(n +2 — 1% “ ten Potenzsumme: 
x Glieder AkS, (mn) =kASG , (mn) =iIni kN). 


x Hierin bilden die ersten n Glieder das alte $7;(n). Die 


folgenden x Glieder können aus den Gliedern von 
Sp(a) Ok +1E +: + (x — 1) gebildet werden, 
indem man die ardall jedes Gliedes um n vermehrt; 
genau das drückt man aus, wenn man symbolisch 2) 
für die Summe der fraglichen « letzten Glieder 


(S(z) + n)k 
schreibt. Man erhält so: 


Sk(n + x) = Sp(n) + (S(z) + n)k 


1) vgl. z.B. N.E. Nörlund, Differenzenrechnung. Berlin 
1924, 8. 17/21. 

2) "Die Potenz ist nach dem binomischen Lehrsatz auszuführen. 
Dann sind die Exponenten von B (bzw. S) als Indices zu schreiben, 


Daher hat die Differenz des Unterschiedes beider 
Funktionen für alle n den Wert 0: 
A XS, (n) -k Sm} = 0. 

Der Unterschied selbst muß demnach eine von n un- 
abhängige Konstante Bj. sein: 

En) = kei...) 
Die Rekursionsformel (9) ist äußerst bequem, um die 
S%;(n) durch ee auseinander zu gewinnen 3): 

S;(n) = ff Seal) da tt Bun... (105 


8) Diesen Weg een Herr I.Paasche, München, in 
einer brieflichen Mitteilung an mich. 
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denn nach Definition (3) der Sg ist 8%(0) =. Zur Be- 
stimmung des Zahlenwertes von B7; kann wiederum (8) 
dienen. 


Kennt man z. B. schon 


nn —)@r—) _ nm n 
ee 
so folgt aus (10) für k=3: 
N 
Ze 25 s 
sn=3 |(E-8+4 dz + Bzn 
0 


n*  n® n? 
ge 


Hier liefert »n = 1 links gemäß (8) S, (1) = 0% = 0, also: 


1 1 1 
Ten ı 
Es folgt 
4 3 2 2 — 1;)2 
B=0, Sm) =-ı-—- ee 


4 2 4 4 


Die Rekursion liefert also zugleich die Potenzsummen 
und die B7. 


Z. angew. Math. Mech, 
Bd. 34 Nr. 3 März 1954 


Es bleibt zu zeigen, daß die Br die Bern oul l IS 
schen Zahlen sind. Dazu verschafft man sich die 
Taylorkoeffizienten von ;,(n) durch Weiter- 
differenzieren von (9) und Einsetzen von n=0. Man 
erhält die Taylor entwicklung 


S,(n) = Br + () Br_1 a %) Ban? +:-- 


+ (K) Bone, =, 


wobei zu den durch (9) definierten B;(k 1) noch 
B, =1 hinzugenommen wurde. Die Gestalt dieser 
Entwieklung legt symbolische Schreibung nabe: 


Sm)=(B+me.....(iD. 


Einsetzen von n = 1 liefert für k > 2 nach (9) und (8) 
sofort die Definitionseigenschaft (l)der Bern o ull r 
schen Zahlen. Einsetzen von (11) in (9) gibt die 
explizite Darstellung (4) der Potenzsummen durch 
Bernoulli sche Polynome. . 

Im übrigen ist (11) unmittelbar interessant: Die Ab- 
leitung der k-ten Potenzsumme ist danach identisch 
mit dem k-ten Bernoullischen Polynom, eine Be- 
ziehung, die man zwar mit Hilfe bekannter Differen- 
tiationseigenschaften der Bernoullischen Poly- 
nome sofort aus (4) ablesen kann, die aber merkwürdiger- 
weise wenig geläufig zu sein scheint. 


Darmstadt. Th. Zech. 


BUCHBESPRECHUNGEN 


Enzyklopädie der math. Wissenschaften 
mit Einschlußihrer Anwendungen. Leipzig 
1953. B. G. Teubner, Verlagsgesellschaft. i 


G. Pickert, Lineare Algebra. Normalformen 
von Matrizen, Bd. Il, Heft 3, Teil I, 6/7. 728. Preis 
geh. 7,50 DM. 

Die gedrängte Wiedergabe algebraischer Themen 
ist wegen der Abstraktheit des Gegenstandes beson- 
ders groß; trotzdem läßt sie sich überwinden. Die 
allgemein orientierende geschichtliche Einleitung ist 
darum zu begrüßen. Sie ist aber nicht hinreichend 
dafür, daß der Eindruck unmöglich wird, die Ab- 
straktion sei um ihrer selbst willen da. Ein die leben- 
dige Entwicklung widerspiegelndes und die Notwen- 
digkeit der modernen Begriffsbildung aufzeigendes 
Bild vermittelt der Artikel leider nicht. Er ist ledig- 
lich eine Schattenprojektion. 


W. Maak, Darstellungstheorie unendlicher 
Gruppen und fastperiodische Funktionen, 
Bd. I1, Heft 7, Teil L,16. 26 S. Preis geh. 3,— DM. 

Der Artikel deckt sich inhaltlich großenteils mit 
des Vf. Lehrbuch über „Fastperiodische Funktionen“. 
Er gibt nicht nur eine gut orientierende Übersicht 
über Begriffsbildung und Eigenschaften, sondern 
unterrichtet auch über die mannigfachen Beziehungen 
zu benachbarten Problemen und über den historischen 
Werdegang. Der Frage der Anwendungen und An- 
wendungsmöglichkeiten ist ein besonderer Abschnitt 
gewidmet. 


Ph. Furtwängler, Allgemeine Theorie deralge- 
braischen Zahlen Bd. I2, Heft8, Teil II, 19. (Über- 
arbeitet von H. Hasse u. W. Jehne, Hamburg) 50 S. 
Preis geh. 6,— DM. 

Der Artikel beweist, daßauch vom heutigen Stand- 
punkt der Forschung aus ein Enzyklopädie-Artikel 
möglich ist, in dem die geschichtliche Entwicklung 
einer mathematischen Theorie mit allem jetzt ein- 
sichtig gewordenen Beziehungsreichtum verarbeitet 


wird. Dadurch entsteht sowohl für den bloß Nach- 
schlagenden wie für den Leser ein plastisches Bild 
und eine klare Vorstellung von der inneren Not- 
wendigkeit der Weiterbildung der Theorie. Die Rolle 
und Funktion der modernen Abstraktion wird be- 
griffen in des Wortes ursprünglicher Bedeutung. Die 
Vorzüge, die dieser Artikel in sich vereinigt, können 
gar nicht überschätzt werden. 


Dresden. 


H. de Wolf Smyth, Atomenergie undihre 
Verwertung im Kriege. (Aus dem Eng- 
lischen übersetzt von Prof. Dr. F. Dessauer 
Freiburg/Schweiz.) 352 Seiten, Basel 1947. Ernst 
Reinhard Verlag A. G. Preis kart. 11,15 DM; geb. 
14,50 DM. 

Bei dem Buch handelt es sich um den offiziellen Be- 
richt über die Anstrengungen, die während des zweiten 
Weltkrieges in den Vereinigten Staaten von Amerika, 
England und Kanada zur Entwicklung der Atom- 
energie für Kriegszwecke gemacht wurden. Dieses 
Stück der Wissenschafts- bzw. Technikgeschichte hat 
durch seine politischen Folgen so allgemeine Be- 
deutung, daß die sehr ausführliche gemeinverständ- 
liche Darstellung der historischen Entwicklung auch 
heute noch Interesse beansprucht. Der Naturwissen- 
schaftler dürfte beim Lesen dieses Berichtes in erster 
Linie gefesselt werden durch die Erkenntnis, zu wel- 
chen ungeahnten Erfolgen die planvoll gelenkte Zu- 
sammenarbeit der Wissenschaftler führen kann und 
welche Möglichkeiten des technischen Fortschritts, 
aber auch des Mißbrauchs, sich dabei bieten. 


Dresden. Recknagel. 
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Berlin: Am 28. November 1953 verstarb nach 
kurzem Leiden unerwartet Dr.-Ing. E. h. Edmund 
Altenkirch im 74. Lebensjahr. 
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